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Sur les Equations aux Dérivées Partielles de la Physique 

Mathématique. 

Par H. Poincaré. 



Quand on envisage les divers problèmes de Calcul Intégral qui se posent 
naturellement lorsqu'on veut approfondir les parties les plus différentes de la 
Physique, il est impossible dé n'être pas frappé des analogies que tous ces prob- 
lèmes présentent entre eux. Qu'il s'agisse de l'électricité statique ou dynamique, 
de la propagation de la chaleur, de l'optique, de l'élasticité, de l'hydrodynamique, 
on est toujours conduit à des équations différentielles de même famille et les 
conditions aux limites, quoique différentes, ne sont pas pourtant sans offrir quel- 
ques resemblances. Nous ne citerons ici que quelques exemples. 

J'imagine d'abord que l'on se propose de trouver la température finale d'un 
corps solide conducteur, homogène et isotrope, lorsque les divers points de la 
surface de ce corps sont maintenus artificiellement à des températures données. 

Ce problème traduit dans le langage analytique s'énonce comme il suit : 

Trouver une fonction V qui dans une portion de l'espace satisfasse â l'équa- 
tion de. Laplace, 

<Z 2 F d?V d?V_ 
A V * dtf + dtf + di — ° ' 

et qui prenne des valeurs données aux divers points de la surface qui limite cet 
espace. 

C'est le problème de Dirichlet. 

Supposons maintenant que l'on cherche quelle est la distribution de l'élec- 
tricité statique à la surface d'un conducteur donné ; nous retrouverons le même 
problème analytique. 

Il s'agit de trouver une fonction Y qui satisfasse à l'équation de Laplace dans 
tout l'espace extérieur au conducteur et qui se réduise à à l'infini et à 1 à la sur- 
face au conducteur. 
28 
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C'est un cas particulier du problème de Dirichlet, mais on connaît un moyen 
(par les fonctions de G-reen) de ramener le cas général à ce cas particulier. 

Les deux problèmes, absolument différents au point de vue physique, sont 
identiques au point de vue analytique. 

D'autres analogies, quoique moins complètes, sont cependant évidentes. 

Nous citerons d'abord le problème suivant ; un liquide est contenu dans un 
vase qu'il remplit complètement ; divers corps solides mobiles sont plongés dans 
ce liquide ; on connaît les mouvements de ces corps et on suppose qu'il y a une 
fonction des vitesses ; on demande quel est le mouvement du liquide. 

C'est là le problème des sphères puisantes de M. Bjerknes (imitation hydro- 
dynamique des phénomènes électriques). 

Au point de vue analytique, il s'agit de trouver une fonction V qui satis- 
fasse à l'équation de Laplace à l'intérieur, d'un certain espace et telle que sur la 

surface qui limite cet espace la dérivée -=- ait des valeurs données. 

dV 
Je rappelle quel est le sens de cette notation -=— dont il sera fait un fréquent 

QJÏt 

usage dans la suite. Soit un élément de surface quelconque ; et a, (3 , y les trois 
cosinus directeurs de la normale à cet élément ; nous posons : 

dV_ dV dV dV 
dn dx " dy ' dz 

Ainsi dans le problème hydrodynamique, nous retrouvons la même équation 
différentielle que dans les problèmes thermique et électrique ; les conditions aux 
limites seules différent. Il en sera encore de même dans le problème de l'induc- 
tion magnétique. 

Supposons un ou plusieurs aimants permanents mis en présence d'un corps 
magnétique parfaitement doux M. Il s'agit de trouver une fonction V (le 
potentiel magnétique) qui satisfait à l'équation de Laplace dans toute la por- 
tion de l'espace qui n'est pas occupée par des aimants permanents et qui est 
assujettie en outre aux conditions suivantes. Aux divers points où il y a 
du magnétisme permanent, AT n'est pas nul, mais peut être regardé comme 
donné. La fonction V est continue dans tout l'espace ; ses dérivées sont con- 
tinues à l'intérieur du corps M et à l'extérieur de ce corps, mais elles sont 

discontinues à la surface du corps M. Dans le voisinage de cette surface, 

dV 

-j— aura donc deux valeurs différentes selon qu'on se placera à l'intérieur 
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ou à l'extérieur du corps M) mais le rapport de ces deux valeurs sera une con- 
stante donnée. 

Ici encore, nous retrouvons la même équation différentielle, avec des con- 
ditions aux limites analogues quoique différentes. 

Voici maintenant des cas où l'équation différentielle est légèrement modifiée. 

Supposons que l'on cherche la loi du refroidissement d'un corps solide isolé 
dans l'espace. Il s'agira de trouver une fonction V satisfaisant à l'équation 

dV 7 ATT 

W = hAV > 

et qui de plus est donnée pour t = 0. Enfin à la surface du corps le rapport de 

V à -^— est donnée. 
cm 

Dans les problèmes d'optique, on a trois fonctions inconnues u, v, w et 

quatre équations : 

d?u , . d 2 v , . cPw , . du dv dw 

Les conditions aux limites varient suivant les problèmes ; mais dans les questions 
de diffraction principalement elles ne sont pas sans analogie avec celles que nous 
avons rencontrées jusqu'ici. 

Ce sont encore les mêmes équations, avec des conditions aux limites ana- 
logues, quoique différentes, que l'on rencontre dans le problème de la viscosité 
des liquides, traité d'après des idées de Navier. Les récents travaux de M. 
Couette ont rappelé l'attention sur cette question, qui était tombée dans un 
injuste oubli, malgré le beau chapitre que Kirchhoff y avait consacré dans sa 
Physique Mathématique. 

La théorie de l'élasticité nous offre des équations plus compliquées, mais 
qui ne différent pas beaucoup des précédentes. 

On a encore trois fonctions inconnues u, v, w, auxquelles j'adjoindrai la 

fonction auxiliaire 6 = -^ — l- -, — I — =— et trois équations dont la première s'écrit: 

dx dy dz 

Au + ù-^ = 0. 

Je n'écris pas les conditions aux limites tout à fait analogues, mutatis rnutandis, 
à celles des problèmes précédents et je passe immédiatement à une question très 



214 PoincarÉ: Sur les Equations aux. Dérivées 

importante en hydrodynamique et qui consiste à trouver les composantes de la 
vitesse en tous les points d'un liquide quand on connaît les composantes du tour- 
billon en tous les points de ce même liquide. Ce problème au point de vue 
analytique s'énonce comme il suit : 

Connaissant trois fonctions a, (3, y, trouver trois fonctions inconnues 
u, v, w, qui satisfont à certaines conditions aux limites et de plus aux équations 

_dw dv „ du dw dv_ du du dw dw . . 

a -~&y~~àz' P—^dJ~~ûc' 'y—d^~~dy , dx + dy + dz~ U ' W 

L'analogie avec les problèmes précédents ne paraît pas d'abord évidente, mais 
elle le devient si on observe que les trois premières équations (1) peuvent (en 
vertu de la quatrième) être remplacées par trois autres dont la première s'écrit 

A _§§_ dy 

dz dy ' 

et dont des autres peuvent s'écrire par symétrie. 

Je pourrais montrer aussi, si je ne craignais de fatiguer l'attention par de 
trop nombreux exemples, que presque toutes les questions, encore mal étudiées, 
relatives à l'induction électrodynamique dans des conducteurs non linéaires se 
ramènent à des problèmes analogues' aux précédents et surtout à la dernière 
question d'hydrodynamique que je viens de mentionner. 

Cette revue rapide des diverses parties de la Physique Mathématique nous 
a convaincus que tous ces problèmes, malgré l'extrême variété des conditions aux 
limites et même des équations différentielles, ont, pour ainsi dire, un certain air 
de famille qu'il est impossible de méconnaître. On doit donc s'attendre à leur 
trouver un très grand-nombre de propriétés communes. 

Malheureusement la première des propriétés communes à tous ces problèmes, 
c'est leur extrême difficulté. Non seulement on ne peut le plus souvent les 
résoudre complètement, mais ce n'est qu'au prix des plus grands efforts qu'on 
peut en démontrer rigoureusement la possibilité. 

Cette démonstration est-elle nécessaire? La plupart des physiciens en 
feraient bon marché. L'expérience ne permettant pas de douter, par exemple, 
de la possibilité de l'équilibre électrique, on ne peut douter, non plus semble-t-il, 
de la possibilité des équations qui expriment cet équilibre. Nous ne saurions 
nous contenter de cette défaite ; l'analyse doit pouvoir se suffire à elle-même et 
d'ailleurs un pareil raisonnement, n'il s'applique peut-être aux problèmes que 
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l'on rencontre directement en Physique, ne saurait s'appliquer de même à une 
foule de problèmes plus simples, qui se posent d'eux mêmes dès qu'on cherche 
à résoudre les premiers. En outre toute démonstration rigoureuse de la possi- 
bilité d'un problème en est toujours une solution ; dans le cas qui nous occupe, 
cette solution sera généralement grossière et tout à fait impropre au calcul 
numérique; cependant elle nous enseignera toujours quelque chose. 

Maintenant si cette démonstration est nécessaire, devons-nous pourtant nous 
astreindre à la même rigueur que dans une question d'analyse pure ? Ce serait 
dans beaucoup de cas un pédantisme bien inutile. Les équations différentielles 
auxquelles obéissent les phénomènes physiques n'ont été souvent établies que par 
des raisonnements peu rigoureux ; on ne les regarde que comme des approxima- 
tion; les résultats expérimentaux, auxquels il s'agit de comparer les consé- 
quences de la théorie, sont eux-mêmes approximatifs. Dans ces conditions, la 
rigueur absolue est de peu de prix, et il semble souvent qu'il n'y a pas lieu de la 
rechercher si on doit la payer de trops d'efforts. 

Mais alors comment reconnaîtra-t-on qu'on raisonnement dont la rigueur 
n'est pas absolue, n'est pas un simple paralogisme ? Quand aura-t-on le droit de 
dire que telle démonstration, insuffisante pour l'Analyse, est assez rigoureuse 
pour la Physique ? La limite est bien difficile à tracer. J'essayerai pourtant de 
le faire ; je m'efforcerai de marquer nettement cette frontière et d'expliquer 
pourquoi en deçà on est encore dans le domaine de la science, et au delà dans 
celui du paralogisme. 

Néanmoins toutes les fois que je le pourrai, je viserai à la rigueur absolue 
et cela pour deux raisons ; en premier lieu, il est toujours dur pour un géomètre 
d'aborder un problème sans le résoudre complètement ; en second lieu, les équa- 
tions que j'étudierai sont susceptibles, non seulement d'applications physiques, 
mais encore d'applications analytiques. C'est sur la possibilité du problème de 
Dirichlet que Pdemann a fondé sa magnifique théorie des fonctions abélienne. 
Depuis d'autres géomètres ont fait d'importantes applications de ce même prin- 
cipe aux parties les plus fondamentales de l'Analyse pure. Est-il encore permis 
de se contenter d'une demi-rigueur ? Et qui nous dit que les autres problèmes 
de la Physique Mathématique ne seront pas un jour, comme l'a déjà été le plus 
simple d'entre eux appelés à jouer en Analyse un rôle considérable? 
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§1. — Problème de Dirichlet. 

Le Problème de Dirichlet énoncé plus haut est toujours possible ; ce prin- 
cipe est connu sous le nom de principe de Dirichlet et la première démonstra- 
tion est due à Riemann. Si une fonction F est assujettie à prendre des valeurs 
données aux divers points d'une certaine surface, limitant un certain volume. où 
le fonction et ses dérivées sont continues, l'intégrale triple : 

'dTV , fdV\* , fdY\ 
.dy . 

ne peut s'annuler ; elle admet donc un certain minimum et il est aisé de vérifier 
que ce minimum correspond au cas où la fonction Y satisfait à l'équation de 
Laplace. Cette démonstration, qui est à peu de chose près elle de Riemann 
n'est pas rigoureuse car elle est soumise à toutes les objections relatives à la con- 
tinuité des fonctions définies par le calcul des variations. 

Aussi un très grand nombre de géomètres se sont-ils préoccupés d'établir 
plus solidement ce principe ; je citerai en première ligne les recherches de M. 
Schwarz (dans le programme de l'Ecole Polytechnique de Zurich 1869 et dans 
les.Monatsberichte de l'Académie de Berlin 1870) bien qu'elles se rapportent 
plus particulièrement au cas de deux variables et ne puissent pas toujours 
s'étendre sans modification au cas qui nous occupe. 

M. Keumann a donné de son côté une méthode générale qui permet de 
résoudre complètement le problème, si la surface où la fonction Y prend des 
valeurs données- est convexe. Il résout donc le problème de la distribution 
électrique dans le cas d'un conducteur convexe. 

La méthode de M. Robin ne s'applique également qu'aux conducteurs con- 
vexes. Toutefois il y a un certain nombre de méthodes plus ou moins com- 
pliquées connues sous le nom de méthodes alternantes et qui permettent d'étendre 
les résultats au cas d'un conducteur de forme quelconque ou de plusieurs conduc- 
teurs isolés. 

Le problème de Dirichlet se ramène à la recherche des fonctions de Green. 

Soit â trouver une fonction Y qui à l'intérieur d'une surface S satisfasse à 
l'équation de Laplace et sur cette surface prenne des valeurs données. 

Supposons qu'on ait trouvé une fonction U qui soit finie, qui satisfasse à 
l'équation A £7"= 0, continue à l'intérieur de S sauf au point (a, b, c) intérieur à 
S où la fonction U sera infinie, de telle fapon que la différence 

1 



U— 



*S(x — af +{y — bf + {z— cf 
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soit finie. Pour achever de définir la fonction U nous l'assujettirons à s'an- 
nuler en tous les points de S. On voit alors que la valeur de Y au point a, b, c 
est égale à l'intégrale 



/ 



Y-j— da 
an 



An 
étendue à tous les éléments da de S. 

Donc quand on saura trouver les fonctions de Green, on saura résoudre le 
problème de Dirichlet. 

D'autre part la recherche des fonctions de Green se ramène, à l'aide de la 
transformation par rayons vecteurs réciproques, au problème de la distribution 
électrostatique à la surface d'un conducteur. 

Ce problème de la distribution électrostatique n'est qu'un cas particulier du 
problème de Dirichlet, et cependant nous voyons que le cas général s'y ramène. 
La méthode alternante de Murphy permet d'ailleurs de ramener le problème de 
la distribution à la surface de plusieurs conducteurs isolés, au cas d'un conduc- 
teur unique. Nous ne nous occuperons donc plus désormais que du cas par- 
ticulier où l'on cherche la distribution électrostatique à la surface d'un seul 
conducteur, puisque le cas général 1 y ramène". 

Que devons-nous penser des méthodes proposées jusqu'ici ? Ce sont à la 
fois des méthodes de démonstration destinées à établir la possibilité du problème 
et des méthodes de calcul destinées à le résoudre effectivement. Comme 
méthodes de démonstration, elles sont assez compliquées, mais elles se com- 
plètent mutuellement de façon à s'appliquer à tous les cas et à satisfaire les juges 
plus sévères au sujet de la rigueur. 

Comme méthodes de calcul, elles ne valent rien ; car personne n'aura 
jamais l'idée de les appliquer ; même les plus simples d'entre elles, celles de 
Neumann ou de Robin, conduisent à des calculs inextricables dès la seconde 
approximation. Tout ce qu'on peut espérer en tirer, sans un labeur par trop 
écrasant, ce sont des inégalités assez grossières auxquelles sera sommise la capacité 
du conducteur. 

Tel est l'état actuel de la question ; voyons maintenant quel est le but que 
je me suis proposé ; ce qui eût été le plus intéressant, c'eût été de remplacer les 
méthodes de calcul actuel par d'autres moins défectueuses. Je n'ai pu le faire, 
je me suis borné à chercher une méthode de démonstration plus simple que celles 
qui ont été proposées jusqu'ici et directement applicable à tous les cas. 
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Je vais commencer par rappeler succinctement les principales propositions 
dont j'aurai à faire usage dans la suite. 

1°. La fonction de Green, U, relative à une sphère S et à un point P s'obtient 
de la façon suivante. Soit R le rayon de la sphère. 

732 

Prenons sur la droite OP une longueur OQ= jfp) la fonction Casera le 
potentiel de deux masses l'une égale à 1 et placée au point P, l'autre égale à 



— V 7T5 e * pl ac ée au point Q. 



OP 

dU 
2°. La valeur de -=— correspondant à un élément quelconque de la sphère 

/S'est en raison inverse du cube de la distance.de cet élément au point P. Elle 

est égale à jp gpa \ 

B JfP 8 ' 

M désignant le centre de gravité de l'élément considéré. 

3°. Si l'on considère une sphère S et un point P intérieur à cette sphère, le 

potentiel d'une masse électrique égale à 1 et placée au point P sera égale à -^^ 

au point M . 

Imaginons ensuite que cette même masse électrique égale à 1 se répartisse 
sur la surface de la sphère S de telle fapon que la densité .sur un élément quel- 
conque de cette sphère soit en raison inverse du cube de la distance de cet 
élément au point P. Je dis que le potentiel de cette masse ainsi distribuée que 

nous appellerons W sera égale à -=j^ en tout point M extérieur à la sphère et 

plus petit que ^p si le point M est intérieur à la sphère. 

En effet considérons une fonction qui soit égale à la fonction de Green, U, 
définie plus haut à l'intérieur de la sphère et à l'extérieur de la sphère. Cette 
fonction satisfait partout à l'équation de Laplace ; elle est continue sauf au point 
P et sur la surface de la sphère. Au point P la fonction devient infinie et sa 

différence avec ^— reste finie ; sur la surface de la sphère, la fonction elle-même 

^ , v iE 2 OP 2 

reste continue, mais sa dérivée -^— subit un saut brusque égal à „ „,„„ . 

dn i o B.MP 3 
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Nous devons en conclure que cette fonction est égale au potentiel de diverses 
masses électriques distribuées comme il suit : 
1°. Une masse égale à 1 au point P. 

jg» OP* 

2°. Une masse de densité — j — w^ps aux divers points de la surface de S. 

Cette seconde masse n'est autre chose que la masse définie plus haut et dont le 
potentiel était W, mais changée de signe. 
On a donc : 

vp- w =° 

à l'extérieur de S et 

à l'intérieur de S. G. Q. p. d. 

4°. Supposons que nous nous proposions de trouver une fonction Y qui satis- 
fasse à l'équation de Laplace à l'intérieur de la sphère S et qui aux différents 
points M de la surface de cette sphère prennent des valeurs données V°. 

La valeur de cette fonction V en un point P intérieur à la sphère sera 

l'intégrale njg g pu yo 






étendue à tous les éléments do de la sphère. 

5°. Soient w et W la plus petite et la plus grande valeur que puisse 
prendre V° ; ce seront aussi, comme on le sait, la plus petite et la plus grande 
valeur que puisse prendre V. 

Si w est positif, il sera de même de V° et de V. 

D'ailleurs comme MP est toujours compris entre B — OP et E + OP; on 

aura Tr ^ B+ OP rV'do 



v< R+OP P 

^(E— opyj < 



(E — OPfJ értB 
et v E— OP rV°d(ù 

^ (E + OPfJ 4rtB ' 

6°. J'arrive au théorème de Harnack dont je ferai un fréquent usage et qui 

peut s'énoncer comme il suit : 

Si une série 

T 7 i+F 8 -f-....+ V n + .... 
29 
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est convergente en un point intérieur à sphère S, que tous ses termes soient 
positifs à l'intérieur de cette sphère et satisfassent tous à l'équation de Laplace, 
cette série sera uniformément convergente à l'intérieur de toute sphère intérieure 
elle-même à S. 

Soit en effet V° la valeur de V t en un point quelconque de la surface de S; 
d'après l'hypothèse V% sera essentiellement positif de sorte qu'on pourra écrire : 

1 ^ (E — OPfJ éTtB ' Vi > (22 + OPfJ 4nB ' 

Si maintenant nous appelons Y[ la valeur de V t au point 0, centre de la sphère, 
on aura *V*da 



*=/ 



AnW 



d'où: v , R(R+OP) ^ B(R-OP) y, 

* (22— opy > Vi > (r+ opy *' 

ou si le point P est intérieur à une sphère s, concentrique à S et de rayon r<C 22, 

22(22 + r)22(22-r) r 
(22 _ r f Vi > Vi > (22 H- rf *' 



Si donc le point P est intérieur à s, chacun des termes positif de la série 

V 1 +V t + ....+ V n +.... 
sera plus petit que le terme de la série convergente 

"FÎ+PÎ + ....+ K +...'• 

multiplié par le facteur constant 

22(22 + r) 
(22 — rf ' 

Si donc la série converge au point 0, elle convergera en tous les points intérieurs 
à S, et la convergence sera uniforme à l'intérieur de s . 

Il suffit d'ailleurs que la série converge en un point P quelconque pour 
qu'elle converge au point 0; on a en effet : 

r<y (22 + opy 
*^ *22(22— OP)' 

Ainsi si la série converge en un point quelconque intérieur à S, elle convergera 
encore en tous les autres points intérieurs à S, et la convergence sera uniforme 
dans toute sphère intérieure à S.. 
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Ce n'est pas tout ; non seulement la série : 

7 X + 7.+ .... + F.+ .... 
converge uniformément à l'intérieur de s ; mais il en est de même des séries 

EL < ^L 4. + ^» + 

dx dx dx ' ' ' ' ' 

*Vi , d?V 2 , d?V n , 

da? "+* M "*" ^ dx % "+" ' 

dMy dxdy 7" * ' ' * dœcfo/ "'"••••» 



La démonstration est la même ; et en effet la fonction F* étant exprimée sous la 
forme d'une intégrale définie, une dérivée quelconque de "F^ s'obtiendra sous la 

même forme par le procédé de la différentiation sous le signe / ; et on en déduira 

comme plus haut deux inégalités auxquelles cette dérivée devra satisfaire. 
On en conclut qu'à l'intérieur de la sphère S, la somme de la série 

Pi + 7, + ....+ 7. + ... . 

est une fonction finie et continue ainsi que toutes ses dérivées et que cette fonc- 
tion satisfait à l'équation de Laplace. 

Le théorème démontré pour une sphère, s'étend aisément à une région R de 
l'espace. 

Si dans la région R les fonctions 

'1) Kg, ... . V n , . . . . 

sont positives et satisfont à l'équation de Laplace. 
Si en un point de cette région la série 

7 X + 7 2 + ....+ F. + .... 

converge, elle convergera dans toute la région et la somme sera une fonction 
finie et continue qui satisfera à l'équation de Laplace. 

Ces préliminaires posés, je puis aborder le problème que j'ai en vue. Il 
s'agit de démontrer la possibilité de l'équilibre électrostatique à la surface d'un 
conducteur isolé. 

Je supposerai qu'en chaque point de la surface de ce conducteur il y a un 
plan tangent déterminé et deux rayons de courbure principaux déterminés. Ces 
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conditions restrictives ne sont pas toutes indispensables. Je préfère pourtant 
me les imposer d'abord et rechercher ensuite, par une courte discussion, quelles 
sont celles dont je puis me débarrasser. 

On peut évidemment" toujours trouver une sphère S telle que le conducteur 
soit contenu tout entier à l'intérieur de cette sphère. 

On peut ensuite trouver un système de sphères en nombre infini 

«i, /b 2 , . . . . , o n , . . . . 

jouissant des deux propriétés suivantes : 1°. Chacune de ces sphères sera tout 
entière extérieure au conducteur. 2°. Tout point de l'espace extérieur au con- 
ducteur appartiens au moins à une des sphères du système. 

Cette proposition paraîtra presque évidente à qui voudra se donner la peine 
d'y réfléchir. 

Nous croyons néanmoins devoir la démontrer en quelques mots. 

Il s'agit de démontrer que l'on peut trouver à l'extérieur du conducteur une 
infinité de points 

Gj, Gg, . . . . , G n 



y m 



tels que pour tout point M extérieur au conducteur, il y ait un point O t plus 
rapproché du point M que de la surface du conducteur. Et en effet s'il en est 
ainsi et que du point 0^ comme centre on décrive une sphère S t ayant pour rayon 
la plus courte distance de O i à la surface du conducteur, la sphère 8 t sera tout 
entière extérieure au conducteur et le point M sera intérieur à la sphère S { . 

Cherchons donc à établir l'existence des points G t . 

Considérons d'abord une région finie B située à une distance finie de la sur- 
face du conducteur et soit è la plus courte distance de la région R à cette surface. 
Cela posé considérons le triple système de plans : 

X_ V3' V-VS' " = V3 

parallèles aux trois plans de coordonnées. Les trois quantités m-^ , ?n 2 , w? 3 sont 
des nombres entiers positifs ou négatifs. Ces plans partageront l'espace en une 
infinité de cubes égaux dont la diagonale est égale à S. 

Ceux de ces cubes qui appartiendront en totalité ou en partie à la région E 
seront en nombre fini. Soient 
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ces cubes. L'intérieur de chacun d'eux j'envisagerai un point appartenant à la 
région B; j'obtiendrai ainsi n points 

Oj, Og, . . . . , C n , 

le point G t appartenant à la fois au cube K t et à la région B. 

Si nous considérons maintenant un point M quelconque de la région B, ce 
point appartiendra à l'un de nos cubes, par exemple au cube K t . La distance 
du point M au point G { sera plus petite que S diagonale du cube ; la distance 
du point G t à la surface du conducteur est elle-même plus grande que 8, puisque 
Gi appartient à B. 

Donc tout point de B est plus rapproché de l'un des n points G lt G % , .... O n 
que celui-ci ne l'est du conducteur. 

Je partage maintenant l'espace extérieur au conducteur en une infinité de 
régions 

.... -ti—n, -tl—n + i, . . . . , £1—2, B_i, Sl , itj, . . . . , M n , .... 

Chacune de ces régions est ainsi représentée par la lettre B affectée d'un 

indice qui peut varier depuis — oo jusqu'à + a». La région B t se compose des 

points dont la plus courte distance à la surface du conducteur est comprise entre 

Q* et Q t+1 . Dans chacune de ces régions je définirai les points G lt G % , . . . , G n 

comme je l'ai dit plus haut. J'obtiendrai ainsi une infinité de points G. Il est 

clair alors que tout point extérieur au conducteur sera toujours plus près de l'un 

de ces points que celui-ci ne l'est lui-même de la surface du conducteur. 

c. Q. p. D. 
L'existence des sphères S t est donc établie. 

Il est clair qu'on pourrait construire ces sphères 8 ( d'une infinité d'autres 
manières et que celle que je viens d'exposer en détail dans le seul but de fixer 
les idées n'a été choisie d'une fapon tout à fait arbitraire. La seule chose essen- 
tielle c'est que les sphères 8 it quoique en nombre infini, puissent être rangées 
en une série linéaire : 

où chacune d'elles ait un indice entier positif parfaitement déterminé, c'est-à-dire 
pour parler le langage de M. Oantor, que l'ensemble des sphères /$• soit de la 
l ère puissance. 

L'ordre dans lequel ces sphères seront rangées dans la série linéaire 

Oi , «g , . . . . , Ù n , . . . . 

est d'ailleurs indifférent. 
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Il s'agit maintenant d'établir qu'il existe une fonction Y finie et continue 
ainsi que toutes ses dérivées à l'extérieur du conducteur et satisfaisant à l'exté- 
rieur du conducteur à l'équation de Laplace ; et de plus que cette fonction tend 
vers quand le point (x , y , z) s'éloigne indéfiniment et vers 1 quand le point 
(a;, y, z) se rapproche indéfiniment de la surface du conducteur. 

Soit le centre et B le rayon de la sphère 2 ; imaginons une certaine quan- 
tité d'électricité positive distribuée uniformément à la surface de cette sphère 

avec une densité — = . J'appelle V le potentiel dû à cette électricité. Nous 



aurons 



F,= l 



à l'intérieur de la sphère 2 et en particulier à l'intérieur du conducteur et 

Y— B 

V ° OM 

lorsque le point M (as, y, z) est extérieur à 2. On a dans tous les cas 

0<F <1 

et V tend vers quand le point M s'éloigne indéfiniment. 

Nous allons faire maintenant une série d'opérations que je vais définir. 

Nous avons vu plus haut que si une masse électrique P se trouve à l'inté- 
rieur d'une sphère S, on peut la remplacer par une masse égale distribuée à la 
surface de la sphère de fapon que la densité en chaque point de cette surface soit 
en rayon inverse du cube de la distance de ce point au point P. Le potentiel 
par rapport à un point extérieur à S n'est pas changé, le potentiel par rapport à 
un point intérieur est diminué. 

On peut appeler cette couche électrique répandue à la surface de S la couche 
équivalente à la masse unique P . 

Cela posé, supposons qu'on remplace toutes les masses électriques qui 
peuvent exister, à l'intérieur d'une sphère S ( par la couche équivalente répandue 
à la surface de cette sphère. Le potentiel en un point extérieur à S t ne changera 
pas, le potentiel en un point intérieur à $ diminuera. Cette opération pourra 
s'appeler: "balayer la sphère $." 

Nous partons de la masse électrique répandue sur 2 et dont le potentiel est 
V . Chacun des points de 2 appartenant à l'une des sphères S i} quelques-unes 
de ces sphères contiendront de l'électricité ; soit $ une de celles qui en contien- 
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nent; "balayons" cette sphère; balayons ensuite la sphère S % si cette sphère 
contient de l'électricité et ainsi de suite. 

Il faut diriger les opérations de fapon que chaque sphère soit balayée une 
infinité de fois. On peut par exemple balayer les sphères dans l'ordre suivant : 

Il est aisé de constater que de cette fapon chaque sphère est balayée une infinité 
de fois. 

Soit ensuite Vi ce que devient le potentiel V après la première opération, 
V 2 ce que devient F x après la seconde opération ; soit enfin V n ce qui devient le 
potentiel après n opérations. 

Supposons que la n e opération consiste à balayer le sphère S k . On aura : 

V n = T^_i à l'extérieur de S & , 
V n <C V n -i à l'intérieur de S k . 

On aura donc dans tous les cas : 

Vn<Vn-l. 

Cette inégalité montre qu'en un point quelconque de l'espace V n est toujours 
décroissant (ou du moins toujours non-croissant) quand on fait croitre l'indice n. 

Il importe de remarquer qu'il n'y a aucun moment et en aucun point de 
masse électrique négative. Au début nous avons sur la sphère 2 une couche 
électrique uniforme et positive. Aucune des opérations subséquentes ne peut 
introduire de masses négatives. En effet le balayage d'une sphère quelconque 
consiste à remplacer les masses électriques positives situées à l'intérieur de cette 
sphère par des couches équivalentes positives répandues à la surface de cette 
sphère. 

On a donc V n > . 

Ainsi en un point quelconque de l'espace V n est toujours positif et décrois- 
sant. Donc quand n croît indéfiniment, V n tend vers une limite finie et déter- 
minée que j'appelle V. J'ai donc démontré l'existence d'une fonction F définie 
en tous les points de l'espace. Il reste à étudier les propriétés de cette fonction. 

Considérons une quelconque des sphères S { . Par hypothèse cette sphère 
sera balayée une infinité de fois. Supposons qu'elle le soit à la af opération, à 

la af, .... à la a® Après chacune de ces opérations, il n'y aura plus 

d'électricité à l'intérieur de S t de sorte qu'on aura : 

AF a = (a* = ai, a 2 , a 3 , .... a n , ... .). 
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Quand a n croît indéfiniment, V an tend vers V de sorte que la série : 

V H + (V at -VJ + (V a -VJ + . . . . +(7 fc -TL.O + . . . . 

converge et a pour somme V. Chacun de ces termes satisfait à l'équation de 
Laplace ; de plus chacun d'eux est négatif (sauf le premier) car on a : 

K n < V^ si comme on le suppose a n >a„_ a . 

Donc en vertu du théorème de Harnack la convergence est uniforme et les séries 
déduites de la précédente par différentiation convergent aussi uniformément. 
Donc à l'intérieur de S t la fonction V est continue aussi que toutes ses dérivées 
et satisfait à l'équation de Laplace. 

Mais par hypothèse tout point de l'espace extérieur au conducteur appar- 
tient au moins à l'une des sphères Si- 

Donc en tout point extérieur au conducteur, V et ses dérivées sont con- 
tinues et l'on a : AF= 0. 

Je dis que F tend vers quand le point M s'éloigne indéfiniment. En effet on a : 

7 >7>0. 

Or V tend vers quand M s'éloigne indéfiniment, donc il en est de même de V. 

Il reste à démontrer que V tend vers 1 quand le point M se rapproche indé- 
finiment de la surface du conducteur. 

Soit donc M un point de la surface du conducteur et suppose que le point 
M se rapproche indéfiniment de M . Par hypothèse il y a au point M un plan 
tangent déterminé et deux rayons de courbure principaux déterminés. On peut 
donc construire une sphère a tangente à la surface du conducteur en M et dont 
le rayon r est assez petit pour que cette sphère soit tout entière contenue à 
l'intérieur du conducteur. 

Soit le centre de cette sphère. La fonction 

M O r 
MG~ MO' 

regardée comme fonction des coordonnées x, y, z du point M satisfait à l'équa- 
tion de Laplace et si réduit à 1 à la surface de la sphère a . 

D'autre part la fonction V n est un potentiel dû à des masses électriques qui 
sont toutes positives. Il en résulte que cette fonction V n peut avoir des maxima, 
mais ne peut avoir de minimum. 
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De plus V n est égal à 1 en tous les points intérieurs au conducteur et en 
particulier sur toute la surface de a . En effet cela est vrai de V , mais, en 
balayant l'une des sphères S i} on ne change pas le potentiel à l'extérieur de cette 
sphère, ni par conséquent à l'intérieur du conducteur qui est tout entier extérieur, 
par hypothèse, à toutes les sphères S t . On a donc à l'intérieur du conducteur 

1= V -V 1 =....= V n =.... 

La différence TT _ -^ __ r 

" n MO 

est un potentiel dû aux masses électriques qui engendrent V n et qui toutes sont 
positives et extérieures à a , et à une masse — r concentrée au point Cet par 
conséquent intérieure à a. Toutes celles de ces masses qui sont extérieures à a 
sont positives. Donc à l'extérieur de a, la fonction U ne peut avoir que des 
maxima et pas de minima. À la surface de a , Î7est nul, car 



n MO 
A l'infini, U est encore nul, car 

V =-?L=o 

n MO 

Donc à l'extérieur de a, Une peut être que positif, en sorte que l'on a: 



V n > 



r 



MO' 
Nous avons donc la double inégalité 

T 

i > y» > -jïq » 

qui à la limite devient : 



= ' =MO=r + MM ' 
Il est donc clair que quand la distance MM tend vers , V tend vers 1 . 

c. Q. F. D. 
La fonction Y satisfait donc bien aux conditions que nous nous étions 
imposées et le ^principe de Dirichlet est établi. 

Voyons maintenant si on peut se débarrasser des conditions restrictives que 
nous nous sommes imposées, à savoir que le plan tangent et les rayons de cour- 
bure principaux sont déterminés. Nous ne nous sommes servis de ces condi- 
30 
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tions que pour établir l'existence de la sphère <?, tangente au conducteur et 
toute entière intérieure à ce conducteur. Cette . sphère <y n'existera évidemment 
que si le plan tangent est déterminé ; mais elle pourra exister encore en un 
point singulier où les rayons de courbure ne varieraient pas suivant les lois 
habituelles ; nous n'aurions alors rien à changer à notre démonstration. 

Nous n'avons donc qu'à examiner les points singuliers pour, lesquels la sphère 
a n'existerait pas. S'il y a sur la surface du conducteur de pareils points singu- 
liers, nous ne pouvons pas encore affirmer que la fonction V tende vers 1 quand 
le point if se rapproche de l'un de ces points singuliers; remarquons toutefois 
que l'existence de ces points singuliers n'empêche pas V de tendre vers 1 quand 
M se rapproche d'un point non singulier de la surface. 

Mais on peut aller plus loin ; supposons que le principe de Dirichlet soit 
démontré pour un certain conducteur O. Supposons ensuite que le conducteur 
donné G présente un point singulier Jf , et que l'on puisse construire un con- 
ducteur G", semblable à G', dont la surface passe par if fr et qui soit tout entier 
intérieur à G. Je dis que quand le point if se rapprochera de M , F tendra 
vers 1. 

En effet le principe de Dirichlet établi pour G' l'est aussi pour G" ; il existe 
donc une fonction V" satisfaisant à l'équation de Laplace à l'extérieur dé G" et 
se réduisant à à l'infini et à 1 à la surface de G". Nous allons faire jouer alors 
au conducteur G" 'le même rôle qu'à la sphère G dans le cas précédemment 
examiné. On verrait comme plus haut que : 

l>Vn>V» 

et par conséquent que : 1 >■ V > V". 

Or quand M tend vers M Q : lim V" = 1 . 

Donc lim V = 1 . c. q. f. d. 

Parmi ces points singuliers les plus intéressants sont les points coniques 
que nous allons étudier de plus près. 

Je suppose qu'une partie de la surface du conducteur G soit une portion de 
cône de révolution. Soit i/ le sommet de ce cône. Je dis que F tendra vers 1 
quand le point M se rapprochera indéfiniment de if . 

Mais pour cela il est nécessaire de démontrer le lemme suivant : 

Soit S une sphère fixe, G un cercle fixe sur cette sphère, P un point fixe 
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en dehors de la sphère. Soit maintenant E un ellipsoïde de révolution variable 
assujetti à restes tangent à la sphère S tout le long du cercle G. Soit V une 
fonction satisfaisant à l'équation A7=0 à l'extérieur de E et se réduisant à 
à l'infini et à 1 à la surface de E. Soit V la valeur de V au point P. 

Je dis que V tendra vers 1 quand le grand axe de E croîtra indéfiniment. 

Désignons en effet par p et Vp 2 — 6 2 les axes de l'ellipsoïde E, par p + h 
et \/(p + 7t) a — b % les axes de l'ellipsoïde E', homofocal à E et passant par P. 
On aura : 

p + 6 

A la limite les ellipsoïdes E et I?' se réduisent à deux paraboloïdes homofocaux 
P x et P{ faciles à construire. 

Alors h tend vers une limite finie qui n'est autre que la distance des sommets 
des deux paraboloïdes ; p — b tend vers une limite finie qui est le demiparamètre 
de P x . On voit aussi que i(p — 6) et L(p + h — b) tendent vers des limites 

finies, que L f 1 H — itt) > ten< ^ vers ° et c l ue L (p + &) crort indéfiniment. Par 

conséquent F tend vers 1 . c. Q. F. d. 

Revenons maintenant à notre conducteur G dont une partie de la surface 
appartiendra à un cône de révolution de sommet M Q . Par le point M , et en 
dehors du cône faisons passer une droite M P et supposons que le point M se 
rapproche de M en suivant cette droite ; je dis que V tendra vers 1 . 

En effet du point M je mène une normale au cône et par le pied de cette 
normale je fais passer un parallèle du cône de révolution. Je construis ensuite 
un ellipsoïde de révolution E qui soit tangent au cône tout le long de ce parallèle. 
Je ferai varier cet ellipsoïde de telle fapon qu'il reste constamment tout entier 
à l'intérieur du conducteur G et son grand axe reste fini quand le point M se 
rapproche indéfiniment de M . Cela est manifestement toujours possible. 

Je construis ensuite un potentiel W qui se réduise à 1 à la surface de E. 
Soit W la valeur de W au point M, on aura : 

l>"F>TFo- 
Construisons maintenant une figure homothétique de la précédente en prenant 
pour centre d'homothétie le point M , le cône de révolution sera son propre 



230 Poinoaeé : Sur les Equations aux Dérivées 

homothétique ; l'homothétique de l'ellipsoïde E sera un ellipsoïde E' tangent au 
cône tout le long d'un parallèle ; l'homothétique du point M sera un point M 1 
de la droite M P. Nous choisirons le rapport d'homothétie de façon que M' soit 
fixe ; alors l'ellipsoïde de révolution E sera tangent au cône tout le long d'un 
parallèle fixe. 

Soit maintenant un potentiel W qui se réduise à 1 à la surface de E'. Soit 
Wô la valeur de W ' au point M' ; on aura : 

W =Wi. 

Lorsque le point if se rapprochera de M , le rapport d'homothétie et par consé- 
quent le grand axe de E' croîtront indéfiniment. Donc d'après le lemme TP ' 

tendra vers 1 . Donc W et par conséquent V tendront aussi vers 1 . 

c. Q. F. D. 

Cela posé considérons maintenant un conducteur G quelconque, et un point 
singulier M de ce conducteur ; je supposerai qu'en ce point i¥ le plan tangent 
est déterminé, ou que ce point M est un point conique ordinaire. Je pourrai 
alors construire un conducteur G" formé par exemple d'une sphère et d'un cône 
de révolution circonscrit, ayant pour sommet M ; je pourrai choisir l'angle du 
sommet du cône et le rayon de la sphère assez petits pour que G" soit tout entier 
intérieur à G. 

Le principe de Dirichlet est démontré pour G" ; G" est intérieur à C; nous 
devons en conclure, comme nous l'avons vu plus haut que V tend vers 1 , quand 
M tend vers M . 

Nous pouvons résumer cette discussion en disant que le principe de Dirichlet 
est établi pour tout conducteur dont la surface est telle que le plan tangent en cliaque 
point est déterminé sauf en un nombre limite de points coniques ordinaires. 

Comme méthode de démonstration, celle que je viens d'exposer ne laisse 
rien à désirer comme méthode de calcul ; elle ne vaut pas mieux que celles qu'on 
a proposées jusqu'ici et même si le conducteur est convexe, elle est inférieure à 
celles de Neumann et de Robin. Mais même à ce point de vue, je ne regrette 
pas de l'avoir publiée ; en effet, ainsi que nous l'avons vu, il ne faut guère 
compter que sur la première approximation si l'on veut pousser plus loin, même 
avec le procédé de Neumann, on serait conduit à des calculs trop rebutants. 
Chaque méthode donnera donc seulement quelques inégalités, il n'est donc pas 
inutile de multiplier les méthodes ; dans chaque cas particulier, un analyste 
habile saura choisir celle qui convient le mieux, ou mieux encore les combiner 
toutes d'une manière convenable. 
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A ce point de vue, la méthode que j'expose ici offrira à cette analyste habile 
d'assez grandes ressources à cause de son élasticité (si j'ose m'exprimer ainsi). 
Le choix des sphères $ , S-. , . . . . S n reste arbitraire dans une très large mesure ; 
d'ailleurs on peut sans rien changer à la méthode remplacer la sphère 2 par 
d'autres surfaces, ou même prendre pour V une fonction potentielle quelconque, 
pourvu qu'elle se réduise à 1 à l'intérieur du conducteur et que les masses élec- 
triques qui l'engendrent soient toutes positives. 

On pourra encore remplacer les sphères S t par d'autres surfaces, pourvu que 
l'on sache construire sur cette surface la couche équivalente à une masse élec- 
trique donnée intérieure à la surface. On conçoit qu'on pourra profiter de toutes 
ces facilités pour adapter le mieux possible la méthode à chaque cas particulier. 

Parmi ces perfectionnements dont la méthode est susceptible, il en est un 
qui me paraît fort important. Nous avons construit les sphères S { de fapon que 
tout point extérieur au conducteur soit intérieur au moins à l'une de ces sphères. 

Imaginons qu'on construisa seulement assez de sphères Si tout entières exté- 
rieures au conducteur, pour que tout point extérieur au conducteur et intérieur h 
2 soit intérieur au moins à l'une de ces sphères. Les sphères S t dont le rayon 
doit être fini dès qu'on est à une distance finie du conducteur, empéiteront 
d'ailleurs sur la région extérieure à 2 . 

Dans ces conditions tout point de l'espace est ou bien intérieur au conduc- 
teur G, ou intérieur à l'une des sphères $ , ou extérieur à 2 . 

Il ne peut être à la fois intérieur à G et intérieur kS t , ou bien intérieur à 
G et extérieur à 2 ; mais il peut être à la fois intérieur à deux ou plusieurs des 
sphères $ , ou intérieur à l'une ou plusieurs de ces sphères et extérieur à 2 . 
Nous avons vu plus haut comment on pouvait remplacer une masse électrique P 
intérieure d'une sphère 2 , c'est-à-dire par une couche électrique répandue à la 
surface de la sphère et dont le potentiel soit égal à celui de la masse P à l'exté- 
rieur de 2 et plus petit à l'intérieur. 

Je dis maintenant qu'on peut également remplacer une masse électrique Q 
extérieure à la sphère 2 par une couche équivalente, c'est-à-dire par une couche 
électrique répandue à la surface de 2 , et dont le potentiel soit égal à celui de la 
masse Q à l'intérieur de 2 et plus petit à l'extérieur. 

En effet soit le centre de 2 et R son rayon ; soient P et Q deux points 
situés sur un même rayon vecteur OP et tels que : 

OP.OQ = R\ 
Supposons que P soit intérieur à 2; Q sera extérieur à 2. 
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Considérons deux masses l'une égale à 1 et placée en P, l'autre égale à 
— \J ~jjt e t placée en Q. Nous avons vu que le potentiel dû à ces deux masses 

était une fonction U, nulle à la surface de 2 , positive à l'intérieur de 2 et néga- 
tive à l'extérieur. 

Soit W le potentiel de la couche équivalente à la masse P répandue à la 
surface de 2. D'après ce que nous avons vu, nous aurons: 

à l'extérieur de 2 : ^p — W-0 -^p <\/ gji mq 

et àl'intérieurdeS: ^p-W= £p-s/% ^>0. 

La lettre If désigne toujours le point de coordonnées courantes ce, y, z. 
Ces égalités montre que l'on a : 

à l'extérieur de 2 : W< V -qj> -ttq , 
à l'intérieur de 2 : W = V jjp jtq • 

Par conséquent la couche équivalente à la masse intérieure 1 située en P est 

aussi équivalente à la masse extérieure v jyp située en Q . 

Il importe de remarquer une différence essentielle entre les deux cas ; nous 
savons que la masse totale de la couche équivalente à la masse 1 située en P est 

égale à 1 et par conséquent plus petite que \J jjp , c'est-à-dire que la masse 

extérieure située en Q. 

Ainsi la couche équivalente à une masse intérieure a une masse totale égale 
à cette masse ; la couche équivalente à une masse extérieure a une masse totale 
plus petite que cette masse. 

A côté opérations dont il a été question plus haut et qui consistent à 
"balayer" l'intérieur des sphères $, nous pourrons introduire une opération 
nouvelle que nous pourrons appeler le balayage de l'extérieur de 2. 

Cette opération consistera à remplacer toutes les masses extérieures à 2 par 
la couche équivalente répandue à la surface de cette sphère. 
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Cela posé on fera une série de balayages successifs en partant de la fonction 
V et en ayant soin de diriger les opérations de telle fapon que l'intérieur de 
chacune des sphères S i} ainsi que l'extérieur de 2, soient balayés une infinité de 
fois. 

Les raisonnements que nous avons faits plus haut sont encore applicables et 
on obtiendra une série de fonctions 

qui convergeront de la fonction cherchée V. 

Seulement l'approximation sera plus rapide, parce que chaque balayage de 
l'extérieur de 2, peut remplacer le balayage d'une infinité de sphères S t qui 
devraient, dans la méthode primitive, remplir l'espace infini extérieur à 2. 

De plus, nous devons observer que dans la méthode primitive, chaque 
balayage laissait subsister la même quantité totale d'électricité dans l'espace ; au 
contraire, dans la méthode nouvelle, cette quantité totale diminue après chaque 
balayage de l'extérieur de 2. 

La méthode nouvelle se prête par conséquent mieux que la première au 
calcul des capacités. 

Nous avons vu que le problème général de Dirichlet se ramène au cas que 
nous venons de traiter. Nous pourrions donc nous dispenser d'étudier directe- 
ment ce cas général. Cependant je vais montrer que la méthode exposée ci-dessus 
y est directement applicable. 

Soit donc G une surface partageant l'espace en deux régions l'une extérieure 
à la surface, l'autre intérieure. Soit U une fonction quelconque, bien déterminée 
en tous les points de G. 

Nous nous proposons de trouver une fonction V: 

1°. Qui est finie et continue ainsi que ses dérivées tant à l'extérieur de G 
qu'à l'intérieur, mais qui peut devenir discontinue sur la surface G elle-même. 

2°. Qui satisfasse à l'équation de Laplace tant à l'extérieur de G qu'à l'inté- 
rieur, mais qui peut cesser d'y satisfaire sur la surface G elle-même. 

3°. Qui tende vers quand le point M de coordonnées courantes x, y, z 
s'éloigne indéfiniment. 

4°. Qui tende vers U quand le point M se rapproche indéfiniment d'un point 
de G, soit par l'intérieur, soit par l'extérieur. 

1 er cas. Supposons qu'on puisse trouver une fonction V finie et continue 
dans tout l'espace ainsi que ses dérivées des deux 1 ers ordres ; qui soit égale à à 
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l'infini et à U à la surface de G et qui soit telle que AV soit constamment 
négatif. 

Une pareille fonction pourra être regardée comme un potentiel dû*à de 

l'électricité répandue dans tout l'espace et dont la densité — -j— - sera partout 

positive. 

Nous n'avons donc encore ici que des masses électriques positives. 

Construirons maintenant une infinité de sphères Si, S s , . . . . , S n , de telle 
façon que chacune de ces sphères soit tout entière extérieure à la surface G et 
que tout point de l'espace, extérieur à G, soit intérieur au moins à l'une des 
sphères St. 

Balayons ensuite, comme il a été dit plus haut, ces diverses sphères S t en 
dirigeant les opérations de telle sorte que chacune d'elles soit balayée une infinité 
de fois. 

Soit V n ce que devient V après la n e opération ; on aura encore, puisque 
toutes les masses électriques sont positives : 

V n + 1 <V n et V n >0 

ce qui montre que quand n croît indéfiniment V n tend vers une limite finie et 
déterminée que j'appelle V. Pour un point intérieur à G on a V n = V . 

On verrait, comme plus haut, qu'à l'intérieur de chacune des sphères S t , et 
par conséquent pour tout point extérieur à G, la fonction V est finie et continue 
ainsi que ses dérivées et satisfait à l'équation de Laplace. 

On verrait également que V tend vers quand le point M s'éloigne indéfini- 
ment. Il me reste à montrer que V tend vers U quand le point M se rapproche 
de la surface G. 

Nous supposerons pour fixer les idées que le point M se rapproche indéfini- 
ment d'un point M de la surface G en restant extérieur à cette surface. 

Nous construirons une sphère a tangente à la surface G en M et de rayon 
assez petit pour être tout entière intérieure à G. 

Le principe de Dirichlet étant démontré pour une sphère, nous pourrons 
construire une fonction W qui à l'extérieur de la sphère satisfasse à l'équation 
ATP = 0, qui se réduise à O à l'infini, et à V Q à la surface de la sphère. 

Quand le point M tendra vers M , la fonction W tendra vers V et par con- 
séquent vers U. 

Comparons les fonctions V n et W. 
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A la surface de c, on a: 

v n =v = w. 

A l'infini, on a : V n = = W. 

Maintenant V n et W sont deux potentiels; le premier est engendré par des 
masses électriques toutes positives et dont quelques unes sont extérieures à a ; 
le second par des masses qui sont toutes intérieures à a . 

Donc à l'extérieur de a la différence V n — W peut avoir des maxima, mais 
pas de rninima ; et comme elle est nulle à la surface de a elle sera toujours posi- 
tive. On a donc : 

T^o > K. > W e * P ar conséquent à la limite V > F> W. 

Quand if tend vers M , V et W tendent tous deux vers U. Donc F tend aussi 
vers U. G. q. p. d. 

Le principe de Dirichlet est ainsi établi pour la région extérieure à (7; on 
le démontrerait absolument de la même manière pour la région intérieure. 

2 e cas. Supposons maintenant qu'on puisse trouver une fonction V Q qui 
soit finie et continue ainsi que ses dérivées des deux premiers ordres, qui se 
réduise à à l'infini et à U à la surface de C. (Nous ne supposons donc plus 
que AV est toujours négatif.) On pourra trouver une fonction V Q satisfaisant à 
ces conditions toutes les fois que la fonction U sera elle-même finie et continue ainsi 
que ses dérivées des deux premiers ordres. 

La fonction F pourra être regardée comme un potentiel engendré par des 

masses électriques répandues dans tout l'espace avec une densité ^r -3 • Seule- 
ment comme A V n'est pas toujours négatif, les masses ne seront pas toutes posi- 
tives. Nous pourrons alors écrire : 

V n = VL TV 

Vq étant le potentiel dû aux masses positives seulement, et vy le potentiel dû 
seulement aux masses négatives changées de signe. 

Soient U' et U" les valeurs de Vq et VI' à la surface de O, on aura : 

U= U' — U". 

VI et V " n'étant engendrés que par des masses positives, on pourra, comme nous 
venons de le voir, construire des fonctions V et V" qui satisferont à l'équation 
31 
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de Laplace à l'extérieur de G et qui se réduiront respectivement à U' et à U" à 
la surface de G . 

La différence V= V — V" 

satisfera à l'équation A F= et se réduira à U à la surface de G. 

Le principe de Dirichlet est donc encore ici établi. 

3 e cas. Il nous reste à examiner le cas où la fonction U n'est plus con- 
tinue ainsi que ses dérivées des deux premiers ordres. 

Bien des méthodes s'offrent à nous pour généraliser les résultats obtenus 
dans les deux premiers cas. 

Je ferai d'abord observer que si la fonction U est elle-même continue et si 
ses dérivées du 1 er ordre ne présentent de discontinuités que le long de certaines 
courbes analytiques tracées sur la surface G, les démonstrations dont j'ai fait 
usage dans les deux premiers cas peuvent se répéter sans qu'on ait rien à y 
changer. 

Passons au cas général ; nous pourrons trouver deux séiies indéfinies de fonc- 
tions 77" 77, 77 • 77' 77/ U' 

qui â la surface de G jouissent des propriétés suivantes: 

1°. Elles sont finies et continues ainsi que leurs dérivées des deux premiers 
ordres. 

2°. Ona: U n+1 >U n , UUi<Ui, W>U n . 

3°. On a : lim U n = lim Ui=- U pour n infini. 

(Cela n'aura lieu que pour les points de G dans le voisinage desquels U est 
continue.) 

Nous pourrons alors construire deux séries de fonctions 

*i> * ii • ' • • i r n> • • • • y y i) y ii • • • • , y n 
telles que A V n = A VI = 

à l'extérieur de G et V n = U n V' n — TJi 

à la surface de G. On aura alors : 

V , * "> V V ,,<? V V ">V 

Nous conclurons de là que V n tend vers une limite finie et déterminée V. 
Le théorème de Harnack montre que l'on a : 

A V- . D'ailleurs V' n > V > V n . 
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Il reste à montrer que F tend vers U quand le point M se rapproche indéfini- 
ment d'un point M de la surface de G. Pour cela, il faut que dans le voisinage 
de M la fonction U soit continue, sans quoi la proposition qu'il s'agit de démon- 
trer serait fausse et n'aurait même pas de sens. 

Nous voulons démontrer que l'on peut prendre M assez voisin de M pour que 

U-e<V<U+e 
quelque petit que soit e. 

Comme au point M , U est continu, U n et XJ' n tendent vers une limite com- 
mune U quand n croît indéfiniment. Nous pourrons donc prendre n assez grand 
Pour que: _ e 

d'où a fortiori 

Nous regardons n comme désormais déterminé, nous pourrons prendre M assez 
voisin de M pour que : 

y n >u n — \ vkul + \. 

On a alors: y> y n>Un _e >n _ £ 

À 

et V<Vl<Ul + -^<U+e. C.Q.F.D. 

§2. — Problème de Fourier. 

Le problème de Fourier a pour objet l'étude du refroidissement d'un corps 
solide rayonnant. J'ai donné de ce problème, dans une note insérée aux Comptes 
Rendus, une solution plus rigoureuse et plus complète que celles qui ont été pro- 
posées jusqu'ici. Bien qu'elle ne soit pas encore entièrement satisfaisante, je 
crois qu'il ne sera pas inutile de la rappeler et de la développer ici ; car elle va 
nous servir de point de départ naturel pour ce qui va suivre. 

Considérons un corps solide homogène et isotrope, isolé dans un milieu indé- 
fini à travers lequel la "chaleur se propage par rayonnement. Soit Fia tempéra- 
ture d'un point du corps ; ce sera une fonction de x , y , z et t; soit la 
température extérieure. On aura, à l'intérieur du corps : 
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et à la surface du corps dV 



dn 



+ hV=0. (2) 



a 3 est un coefficient constant qui dépend de la conductibilité du corps et de sa 
chaleur spécifique. Quant à h c'est un coefficient positif et constant qui dépend 
du pouvoir émissif du corps. 

Le premier point est d'établir l'existence d'une infinité de fonctions auxili- 
aires 77 7T 77 

C/j, u z , . . . . , u n , . . . . 

ne dépendant que de a;, y, a et satisfaisant aux équations suivantes : 
On aura à l'intérieur du corps : 

AU n + k n U n =0 



et à la surface : dU. 



^ + hU n = 0. 



dn 
Les quantités \, \, . . . . , k n , . . . . 

sont des coefficients constants que je supposerai rangés par ordre de grandeur 
croisante et que je déterminerai plus. complètement dans la suite. 

Enfin pour achever de définir la fonction U n , nous ajouterons que l'intégrale 



§TP n dx 



étendue à tous les éléments de volume dt du corps solide doit être égale à 1 . 
Nous allons pour démontrer l'existence des fonctions U n employer une démon- 
stration analogue à celle par laquelle Riemann établit le principe de Dirichlet. 
Soit F une fonction quelconque et posons : 

A =J % F z dv , 

*=*/^+/[(g>(f )'+(£)>• 

L'intégrale A ainsi que la seconde des intégrales de l'expression B sont étendues 

à tous les éléments dt du volume du solide et l'intégrale / F 2 da à tous les éléments 

da de sa surface. 

Supposons que la fonction F soit assujettie à la condition 

A—l. 
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Les deux intégrales de l'expression B ont tous leurs éléments positifs. B ne 
peut donc devenir négatif. B ne peut non plus devenir nul ; en effet il ne pour- 
rait s'annuler que si tous les éléments des deux intégrales étaient nuls à la fois, 
c'est-à-dire si l'on avait : F= 

à la surface du corps et 

dF_dF_dF_ 

dx dy dz 

à l'intérieur du corps. Il faudrait donc que F fût encore nul à l'intérieur du 
corps, ce qui est impossible à cause de la condition : 

A=fF 2 dr=l. 

B admettra donc un minimum absolu. Soit U x la valeur de F qui correspond à 
ce minimum. Le calcul des variations nous donne : 

■^-èA=fu x SU x dr = 0, 
1 s» j, Cttïttj , rfdUi dhU x dU x dèUx . dJJ x dhTJ x \ ., 

Or le théorème de Green nous donne : 

p/dU x dhU x , dU x dèU x , dU x dSU x \ 7 pdVj. , TT 7 /• TTÏXT . , 

de sorte qu'il vient : 

-1 BB =f(h TJ X +^) è U x do —fA U x è U x ck = . 

SB doit s'annuler toutes les fois que SA s'annule. On doit donc d'après l'une 
des régies du calcul des variations, pouvoir trouver une constante \ telle que 

SB — \SA 
soit nulle quel que soit $U X . On doit donc avoir identiquement: 

f^è + hU \) s U ^° —f( A ^ + ^ u i) è U ^ = ° 

ce qui exige que tous les éléments des deux intégrales soient nuls ou que l'on 
ait à la surface du corps 
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et à l'intérieur : AU 1 + Jc 1 U 1 = 0. 

On a d'ailleurs par hypothèse : 

A=J % U?dv=l. 

L'existence de la fonction TJi est donc démontrée. 
On trouve : 

,=,/^ + /[(^)' + (f)- + (f)> 

ou en vertu du théorème de Green : 

B=fu 1 (- 1 + h Ci) do —f TJA U^dr 

ou en vertu des équations qui défininent JTj.: 

B = h x futdk = h 1 . 

T> 

Ainsi \ n'est autre chose que la valeur du rapport — =- pour F= U x . Comme 

ce rapport atteint son minimum pour F— Ui, nous devons conclure que \ est 

le minimum du rapport —j . 

Prenons pour F une valeur quelconque, nous obtiendrons une certaine 

valeur de -j- qui sera plus grande que J^. C'est donc un moyen de trouver une 

limite supérieure de \. 

Faisons par exemple F= 1 . 

vien ' J. = / ck= volume du corps solide, 

B = h J da — hx surface du corps solide. 

Le rapport -£- est donc toujours plus petit que le rapport de la surface, du solide 

à son volume. 

Cherchons encore une autre inégalité. 

Appelons W le volume du corps solide et S sa surface. 

Prenons pour origine des coordonnées le centre de gravité du volume du 
corps, 
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Appelons /l'intégrale Çrfdt 

c'est-â-dire le moment d'inertie du volume par rapport au plan des yz . 

Appelons H l'intégrale /» 

/ çêdiù 

c'est-à-dire le moment d'inertie de la surface par rapport au plan des yz. Soit 
a" la distance du centre de gravité de la surface au plan des yz. Posons : 

J/"= x S=J xda. 
Faisons maintenant : F— ax + 1 , 

a étant une indéterminée ; il viendra : 

A- f(ax + lfck = a?I+W, 

B = h f(ax + lfda + Ca?ck = a 2 (hH+ W) + 2aM h + Sh 

d'où: ^ a 2 (hH+ W) -\- 2aMh + Sh 

kl< a?I+V 

Il faut maintenant choisir a de telle sorte que le second membre de cette inégalité 
soit minimum. Ce second membre admet un maximum et un minimum qui sont 
les racines de l'équation en % : 

M*W =(hH+W- M)(Sh — %W) 

ou mv— a, ( ish + tp 2 + WEh) + sm — m 2 + wsh = o . 

On a donc : . 



z. ^ ISh + W* + WHh */(ISh + WHh + W 2 ) 2 — SHh? + JSPh* — WSh 
hl< 2IW 2IW 

cette limite est inférieure à la précédente. 

Occupons nous maintenant de démontrer l'existence de la fonction U 2 ; soit 
une fonction quelconque F assujettie aux deux conditions suivantes : 

A=l, G=fFU 1 (h=0. 

On pourra choisir cette fonction de façon que B soit minimum ; soit U 2 la 
fonction qu'il faut choisir pour F afin de rendre B minimum. On devra avoir : 

-1 SB =f(hU % + ^ 2 ) hU t da —fAUjUJ* = 
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toutes les fois que : 

±-SA=fu ss 8U,dr=0 W=fuitUiïx = 0. 

On pourra donc trouver deux constantes /L x et h % telles que : 

hB — h£A—Z\hC=0 

quelque soit S U % . On a donc identiquement : 

de sorte qu'on aura à la surface du corps 

et à l'intérieur : 

AU 2 + h 2 U z + k 1 U 1 = 0. 

Calculons \ et h^. Nous trouvons d'abord : 

f UA TJ % fa + h % f U x TJjk + \f TJ&e = . 

0r fUid* = 1 G =fu x U % dn = 0, 

onauradonc: ^ + f UA^d* = . 

Or le théorème de Green donne : 

/ UA Ufr -f UA U,ck =f U x ^ c?o =/ U % ^ do. 

Mais dU % _ , _. dTJ x _ , „ 

— --fiU 2 , -fa- — hU x . 

Il reste donc simplement : 

f U X A U 2 d* =f UA Uj.dr . 

Mais AU 1 — — h 1 U 1 . 

Donc a* =f UA U % dk - — \f Ui TJ&e - . 

Donc %± est nul. 

Il vient donc ; A U 2 + &2 U s = 

d ' où : f UA & s d* + hf Ui&t = . 
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0r A=fui<fa=l. 

Donc: 

^=-/^^=-/^f^ + /[(f)" + (f)' + (f)> 

ou: 

h = hfuido +f$ (—^d* = B. 

Ainsi hi n'est autre chose que la valeur de B qui correspond à F= U%\ \ était 
la valeur de B qui correspond à F= TJ X et au minimum de B, 
Par conséquent on a : 

et d'ailleurs : dU-t , , _. rt . N . „ , 

— — 1- h tZj = a la surface du corps, 

A U z + &2 Z7 2 = à l'intérieur, 

L'existence de la fonction U z est donc démontrée. 
Soit maintenant U 3 une fonction telle que 

A=fui<fo=l, C=fu 3 U 1 dr = 0, D=fu 3 U 2 ck = 

et choisie d'ailleurs de telle sorte que B soit aussi petit que possible. 

On devra avoir : 8B = 

toutes les fois que 

54 = 0, $O=0, W = 0. 

On pourra donc trouver trois constantes % x , Jl a et h 3 telles que l'on ait identique- 
ment : hB - IcthA — 2\5 G — 27^D = . 

Un raisonnement analogue à celui qui précède montrerait que l'on doit 
avoir à la surface du corps d U s , „ 

et à l'intérieur A *7 3 + Te 3 U 3 + \ Ux + /l 2 C7" 2 = . 

On démontrerait ensuite comme on l'a fait plus haut que ^ et ^ sont nuls et que 
Te 3 est la valeur de B qui correspond à F-= U 3 . 
32 
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D'après la définition de U 3 , h 3 est donc la plus petite valeur que puisse 
prendre l'expression B quand la fonction F est assujettie aux conditions : 

fj^ék = 1 ; J*F TT x dfg = ; J *FU % ck = . 

D'autre part le % était la plus petite valeur que pouvait prendre B quand F était 
assujettie aux deux premières de ces conditions seulement. Donc : 

^8 ^ "% ' 

La fonction U 3 est ainsi définie par les conditions : 

™± + hU 9 =0, AU, + k,U, = 0, 

J Uidr =1; fu z Ujék = t /V 3 TTjk — . 

La loi est manifeste ; il est inutile de pousser plus loin ce raisonnement. 
On voit que l'on a démontré l'existence d'une infinité de fonctions : 

U\i C'a, .... , Up, . . . . 
telles que l'on a à la surface du corps 

et à l'intérieur : A U p + \ U p = . 

Les coefficients h p sont des constantes positives et telles que : 

Enfin on a: fu,U&=0 P onvp>q 

et fu*fa=sl. 

Ce raisonnement est sujet aux mêmes objections que celui par lequel Riemann 
établit le principe de Dirichlet. Nous nous en contenterons toutefois pour le 
moment, nous chercherons dans la suite à le rendre plus rigoureux. 

Ces fonctions U p ont été entièrement construites par Lamé dans certains cas 
particuliers, par exemple dans celui de la sphère et celui du parallélipipède. 

Dans celui du parallélipipède rectangle dont les trois dimensions parallèles 
aux trois axes sont 2a, 26 et 2c, l'expression des fonctions TJ V et des coefficients 
h p est particulièrement simple. 
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Nous devons avoir en effet à l'intérieur du corps : 

AU+hU=0 
et en outre : 

dU , , -r-r dU -, t~t 

pour œ = a -T7 + ftt/=0, pourœ= — a -5 hu = 0, 

poury = & -j \-hU=0, pour 2/= — h -5 hU—0, 

pour 3 = c -z — \-hZf=0, pour z = — c -= hU=0. 

Posons : 

£T= sin (ylifc + /u x ) sin (A-jœ + /u 2 ) sin (X 3 x-\- fa). 

Les constantes Xj et fa nous seront données par les deux équations : 

Ai cos (/^a + /Ki) -f- h sin (^a + fa) = , 
yli cos {\a — fa) +. A sin (A x a — ^x) = , 

d'où : tg (\a + fa) = tg (X x a — fa) , 

X7t 

x étant un entier ; il suffira de prendre x = d'où : 

sin (%■& + fa) = sin \x 
et x = 1 d'où : sin (Ax« + jtf x ) = cos X x a3. 

On aura alors pour 

^ = ou -z- Xx + h tg (X x a + /«i) = . 
De même fa, fa, \ et X s seront données par les équations: 

fJs=Oou-^ % z + A tg (V + fi 8 ) = 0. 

Enfin on a : h — 1\ + 7\ + Aj. 

Considérons en particulier le cas de & = ; ce qui correspond au cas où la 
surface du corps est imperméable à la chaleur ; il vient 

tg (\a + fa) = 00 
011 * Axa -f- fa = — + xrt x étant entier : 



246 Poincaré : Sur les Equations aux Dérivées 

ou à cause de ^ t = ou — . 

À 

\ = -î^— m 1 étant entier. 
On a alors : , _ / m? ml m| \ 7t 8 

* _ w " + " _ ^ + ~^ r ; t 

où wj x , ffi 2 , rn % sont des entiers. On aura donc: 

& 1 = 0, ^=4^- 

si 2a est la plus grande des trois dimensions du parallélipipède, c'est-à-dire si 

On trouve ainsi h x = , Ui= const. 

& 2 = j-g- , Z7 2 = const. sin 7^x . 

Supposons maintenant Ji = oo , ce qui correspond au cas où la surface du corps est 
maintenue artificiellement à la température ; on a alors : 

tg (\a + (i 1 ) = 

d'où XjO, + fa = xn, x étant entier, 



ou puisque ^ = ou -^ , 



\a = tn 1 - ir , m 1 étant entier. 

À ■ 



Nous trouvons donc encore : 

A_ w " t ""p - + "^;~T' 

m x , wi 2 et »jg étant entiers. 

Mais toutes les solutions ne sont pas acceptables. On doit avoir en effet 

à la surface du corps ; d'où : 

sin (a^a + fr) = sin (/l 2 & + fi z ) = sin (/t s c + ^ 3 ) = , 

ce qui exige que vn x , îw 2 et rn 3 soient au moins égaux à 1 ; il viendra donc : 
Laissons maintenant de côté le cas du parallélipipède rectangle et revenons au 
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cas général; il est clair que, pour un même corps, \, Tc % h n sont des fonc- 
tions du pouvoir émissif h. Je disque ces fonctions sont toujours croissantes. 
Soient en effet h' et h" deux valeurs du pouvoir émissif; soient pour un corps 
donné, k' n et h'J les valeurs correspondantes de h n , et TJ' n et JJH les fonctions TJ n 
correspondantes. On aura : 
à la surface du corps 

et à l'intérieur du corps : 

aui + ku: = aui' + kui' = o. (4) 

Le théorème de Green nous donne : 

fŒ W-^ Ul)d« = f(Ul>AU>- UlAU>:)d«. 

dW dU" 
Dans cette identité remplaçons AC n ', AU", j 1 ) , — par leurs valeurs tirées 

des équations (2) et (4), il viendra : 

(h' - h»)fuiui'd a = (K - wfuiui'â* . 

Supposons que h' et h" diffère très peu de telle sorte que 

TJ' n = TJ'J à des infininients petits près, il vient : 

dhf UL'da = dh n f Ui'dr . (5) 

Dans l'équation (5), les deux intégrales sont essentiellement positives ; il 
reste donc: jt 

dh ^ U ' 

ce qui signifie que h n est une fonction croissante de h. o. Q. F. D. 

Pour h = ; on a évidemment : 

\ = , U x = const. = — r^ ( W étant le volume du corps) . 

En effet, ces quantités satisfont aux équations 

AU 1 = AU 1 + 7c 1 U 1 =0, ^ = 0. 

an 
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L'équation (5) devient alors 

dh ~ W 

S désignant la surface du corps et W son volume. Donc si h est très petit, Je\ 
est sensiblement égal à g 

Nous venons de trouver : 
Nous avons d'autre part : 

,.=,/^ + /[(f)' + (f)V(f) ! ]* 

d '° Ù: K>hfuid<*, fjJld^<~, 

^ K_ dh 
h dk„ 
et dk n j. dh 

h 

Cette inégalité montre que le rapport ■—■ va toujours en décroissant. 

/fc 

Il importe de remarquer, avant d'aller plus loin que quand h est positif, Jc n 
est essentiellement positif. Sans doute, cela résulte de la fapon dont les fonc- 
tions U n ont été définies plus haut ; mais on pourrait imaginer qu'il existe des 
fonctions U autres que celles que nous venons de définir et pour lesquelles on 

aurait : dU P 

^J r j l U=0, 'AU+JeU=0, ! V*ék=l, Je<0. 
dn *J 

Je dis que cela est impossible et il me suffit pour l'établir de montrer que l'on a : 

,=,/^ + /[(f) + Q' + (f)>, 

ce qui se démontrerait par le même calcul que plus haut. 

Si au contraire h était négatif, Je n pourrait aussi devenir négatif. 

Il peut arriver quand on fait varier h, que deux des nombres h n et Je n+1 
viennent à se confondre. Qu'arrivera-t-il alors en général ? 

Soient Je 1 et Je" deux valeurs de Je, U' et U" les fonctions U correspondantes. 
Imaginons que Je 1 , Je", U' et U" soient des fonctions continues de h. Quand 
h<C.h , on aura par exemple Je! <C Je" ; pour h = h , on aura Je 1 = Je" ; pour h > 7i , 
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on aura Té > k". D'après ce que nous venons de démontrer, Té et k" sont deux 
fonctions croissantes de h. 

Maintenant, supposons qu'il y ait n — 1 nombres k inférieurs à la fois à Jsf et 
à y. Nous devrons, d'après nos conventions, appeler k n la plus petite et k n + 1 
la plus grande des deux quantités k 1 et k". Nous aurons donc : 

y =.k n W = k n + i pour h<^h 
et k" — k n k! — k n + 1 pourA>^ . 

Les fonctions k! et k" étant continues et croissantes, les fonctions k n et k n + 1 
définies comme nous venons de le faire seront aussi continues et croissantes. 

En résumé, dans tous les cas possibles, k n est une fonction croissante de h; 
k n atteint donc sa plus petite valeur pour h = . 

Nous allons donc étudier la valeur de k n pour h = . 

Décomposons le volume de notre corps solide, d'une manière quelconque, 
en p volumes partiels. Considérons chacun . de ces volumes partiels comme un 
corps solide de même conductibilité que le solide donné et dont la surface est 
imperméable à la chaleur. Soient : 

^11 1 Ul& >••••> ^ïn 1 ' • • • 1 
Uzi I ^22 )••••) iJsn >••••) 

les fonctions U relatives à ces^j solides. 

fooient : k^ , k^ ,...., k± n ,...., 

%l I "%2 )••••) K%n >•••■) 
"îpl > "p% >••••> ™pn )••••> 

les nombres k correspondants. 

Comme pour chacun des p solides partiels onaA = 0, on aura : 

"ïi — «ki — • • • • —- *î>i — 

et les^> fonctions Î7 U , JZ^, . . . . , U vl seront des constantes. 

Je conserverai la notation JJ n et k n pour les fonctions relatives au solide 
total. Posons maintenant : 

V= a 1 U 1 + a< i U' i + + a n U n , 

les a étant des coefficients indéterminés. 
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Nous avons à la surface du corps, puisque h est supposé nul : 

— -0 
dn 

et par conséquent en vertu du théorème de Green : 

/[(4D'+(f)'+(f)>=-/>^- 

Le second membre peut s'écrire (en vertu de l'équation A Ut + h U t = 0) : 
J (c^CZi +a z U % + ....+ a n U n %<iikiUi + a 2 hU % +.... + aJc n U n ) dt 

ou encore f en vertu des relations / Uldr = 1 , / #î ^««fo = j : 

7^a| + ^a| + + &„a® . 

D'autre part on a : 

fVâ* =f(a l U l + 0,17, + .... + a^tfo = a* + «| + + o», 

d'où 



»î + «a + + «» 



/ 



T 2 ^ 



<&»• (6) 



Nous pouvons disposer de nos n coefficients arbitraires a x , o^, . . . . , a n de fapon 
à satisfaire à, n — 1 conditions. Voici comment nous choisirons ces conditions. 
Nous annulerons d'abord les \ intégrales 



fvU u ck, fvu n ck fvU Ul ck 

étendues au premier solide partiel. 

Nous annulerons ensuite les %% intégrales 

fvUadv fvU^oW 

étendues au second solide partiel. 

Nous annulerons enfin les \ intégrales 

fvu pl ck, ..... fvu px ck 

étendues au dernier solide partiel. 
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Nous devrons avoir d'ailleurs : 

*i + a» + -\-% v — n—\. (7) 

Rappelons quelle est la définition de &i. Xi+1 . Ce coefficient est le minimum du 
rapport des deux intégrales : 

quand la fonction arbitraire Y est assujettie aux /^ conditions : 

fvU n ék —fw^ft = =fvU lx ck = 0. 

Toutes les intégrales sont supposées étendues au premier solide partiel. Donc 
le rapport des deux intégrales B et A étendues à ce premier solide est plus 
grand que \. Kl + 1 . 

On verrait de même que le rapport de ces deux intégrales étendues au 
second solide est plus grand que & 2 .A a +i> etc.; que le rapport de ces deux inté- 
grales étendues au p e solide partiel est plus grand que k p ,^ +1 . 

Donc le rapport des deux intégrales B et A étendues au solide total, sera 
plus grand que la plus petite des p quantités : 

Mais l'inégalité (6) exprime que ce même rapport est plus petit que h n . 

Donc h n est plus grand que la plus petite des quantités (8) pourvu que la 
relation (7) soit satisfaite. 

Ce résultat peut s'énoncer comme il suit : 

Rangeons les quantités 

«ii > ««>••••> "in > • • • • 
%1 ) ! "* > • • • • 



"Î>1 ) > "'put 



par ordre de grandeur croissante. (Il va sans dire que dans la série ainsi obtenue, 
plusieurs termes pourront être égaux ; c'est ainsi que les p premiers termes qui 

sont & u , h n , , h pl sont tous égaux â 0.) 

Chacun des termes de la série ainsi obtenue sera plus petit que le terme 
correspondant de la série : 

Kl, A/ 2 , . . . . , k n , . . . . 
33 
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Cela posé, supposons que notre solide soit un polyèdre limité par des faces 
qui soient toutes parallèles à l'un des trois plans de coordonnées. 

Nous pourrons décomposer notre polyèdre en n — 1 parallélipipèdes rect- 
angles. Soient j. % z. (q\ 

les valeurs de &j correspondant à ces n — 1 parallélipipèdes ; d'après ce que nous 
venons de voir, la valeur de k n correspondant au polyèdre total, sera plus grande 
que la plus petite des quantités (9). 

Si les trois dimensions d'un parallélipipède sont 2a, 2b, 2c de telle sorte que 
a^-by-c; nous avons vu qu'on a, pour ce parallélipipède : 

Si donc aucune des trois dimensions d'aucun de nos n — 1 parallélipipèdes par- 
tiels n'excède a, on aura pour le polyèdre total 

Quand n croît indéfiniment, on peut faire tendre a vers ; donc k n croît indéfini- 
ment. 

Cela est vrai pour un polyèdre formé comme je viens de le dire ; et comme 
on peut construire de pareils polyèdres qui diffèrent aussi peu qu'on le veut d'un 
solide quelconque, on pourrait dire que cela doit être vrai aussi d'un solide quel- 
conque. Mais un semblable raisonnement ne saurait nous contenter. 

Pour démontrer le théorème pour un solide quelconque, je dois d'abord 
chercher une limite inférieure de h % pour un solide convexe quelconque. Par 
solide convexe, j'entends un corps tel que le segment de droite MM, qui joint 
deux points M et M' intérieurs au corps, soit toujours tout entier intérieur au 
corps ; ou ce qui revient au même tel qu'aucune droite ne rencontre la surface 
du corps en plus de deux points. 

Rappelons d'abord la définition de h^; h % est le minimum du rapport 

quand la fonction V et assujettie à la condition : 

fvck—0. (10) 
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On peut transformer la définition dé fapon à faire disparaître cette dernière 
condition. 

Envisageons l'intégrale 

f{V— Vf chat'. 

Dans cette intégrale, dt et dt 1 sont deux éléments de volume quelconques 
du solide donné ; x, y, z, et x 1 , y', z' sont les coordonnées du centre de gravité de 
chacun de ces deux éléments ; F et F sont les valeurs de la fonction V aux 
points x, y, z et x 1 , y', z'; enfin l'intégrale est étendue à toutes les combinaisons 
de deux éléments de volume dt et dt'. (Chaque combinaison sera répétée deux 
fois de la façon suivante ; une première fois le premier élément jouera le rôle de 
dt et le second celui de dt' et la seconde fois ce sera le contraire.) 

On trouve aisément en se servant de la relation (10) et en appelant Wle 
volume total du corps 

f{V— Vf dtdt' = 2 wf V*ck 

de sorte que h% sera le minimum de l'expression 

nm fc/dV\* , /dV\* , /dV\* ll 

2F i[(^-)+w)+(-^-)]^ (11) 

f(V— Vf dtdt'. 

Mais alors la condition (10) devient inutile ; si en effet on ajoute à V une con- 
stante quelconque, la condition (10) cesse d'être satisfaite et l'expression (11) ne 
change pas. 

Nous pouvons donc dire finalement que h% est le minimum de l'expression 
(11), la fonction V étant tout à fait quelconque et n'étant assujettie à aucune 
condition. 

C'est de la transformation de cette expression (11) que nous allons main- 
tenant nous occuper. 

Posons à cet effet : 

x = £ + p cos 4> sin 6, x 1 = % + p' cos ç> sin 6, 

y=V + p sin <£ sin0, y' = y + p' sin q> sin 6, (12) 

z = pcos0 , z' = p'cos6. 

La théorie de la transformation des intégrales multiples nous donne : 
y*( V— Vf dtdt' =f( V— Vfdxdydz dx'dy'dz' 

=f(V— Vf(p — p'f sin cos dd£dyidpdp'd6d<l>. 
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Il va sans dire que, bien que je n'emploie qu'un seul signe / , il s'agit ici d'inté- 
grales sextuples. 

Parlons maintenant des limites d'intégration ; supposons que £ , vt , d et q> 
soient considérés un instant comme des constantes ; quand on fera varier p, le 
point x, y, z décrira une droite ; comme le corps est convexe, cette droite rencon- 
trera la surface du corps en deux points. Soient p et p x les valeurs de p qui 
correspondent à ces deux points d'intersection. 

Pour obtenir alors tous les points (x, y, z), (x!, y', z') intérieurs au corps il 

faudra faire varier p et p' de p à p 1( de à ~ , 4> de à 2n, et donner à £ et 

à v> toutes les valeurs telles que p et p x soient réels. 

Cela posé, cherchons à transformer le numérateur de l'expression (11). 
Nous avons d'abord, en vertu des équations (12): 

dV . . dV , . . . dV , a dV 

—— = cos à) sin -= 1- sin d> sm 6 — = \- cos v -=— 

dp T dx T dy dz 

d'où: 



//(■ 



dV\% 

-y— ) sin ddBdty 



= / J f cos $ sin 6 -5 p- sin $ sm 6 -= — \- cos 6 ~j- \ sm 0d6dq> 



Calculons cette intégrale double en intégrant entre les limites et 2n pour $ et 
entre les limites et -^- pour 6. Le coefficient de ( -=— J sera : 



ffcos*<l> sm 3 6ddd<l> = Ç . 



Il est aisé de voir que le coefficient de (~r~) a la même valeur. 



2n 
3 



Le coefficient de (-5—) sera : 

C Ccos* 6 sin 6 dddxp=: 

„ . . , a dV dV 
Le coefficient de 2 — ; — 5— sera : 
dx dy 

j j cos $ sin 4> sin 3 d0 eZ<?> = . 
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Ceux de 

n dV dV . , a dV dV 

2 -r — et de 2 -^ — seront : 

dx dz dy dz 

J j cos <£ sin 3 cos 6 de d$ =J J sin q> sin 2 6 cos 6 d$ dty = . 
Il reste donc simplement : 

//(£)*.**=S[(40 + (S) , + (£)']- 

Soit if une fonction quelconque de ce, y etjs; la théorie de la transformation des 
intégrales multiples donne : 

/ Mck =J Mdxdydz =J M cos d£ dy dp . 
Il vient donc : 

Posons maintenant pour abréger : 

l'expression (11) deviendra: 

o-rrr I Bsind cos 6 dgdyidddty 
71 J A sin cos OdÇdvidd dty 

Les quantités sous les deux signes / sont essentiellement positives puisque 
varie entre et -^-. Pour trouver une limite inférieure de l'expression (11), 

il nous suffira de connaître une limite inférieure du rapport — -p . C'est de quoi 

nous allons maintenant nous occuper. 

Si la fonction V est choisie de telle façon que A = 1 , il est clair que l'inté- 
grale B ne pourra pas s'annuler ; elle admettra donc un certain minimum. Cher- 
chons à déterminer ce minimum par le calcul des variations. 
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Nous trouvons : 

1 .„ p*dV dôV _ 

L U =ff( V— V')(8V- 3F)(p — p'fdpdp' = 0. 
Transformons ces deux expressions, nous trouvons, par l'intégration par parties, 

>=[f^:-rf>4=o. 

D'autre part, si dans l'intégrale double 

J=ff(v-v)W( 9 - P >yd P d p i 

on permute p et p', l'intégrale ne doit pas changer il vient ainsi : 

J =ff( y, -v)tr( f r-pyapd f r, 
d ' où ±u = *j 

Gt J=ff(V-V')(p-p'fSVd p d p '=0. 

Cela peut s'écrire : 

J=f 1 mVdpz=0 

en posant: 5= r£\p - pQ^ _£> (p-pT ^. 

Posons donc pour abréger : 

L <¥ = Pi" Po; £ pW = £l_Po . ^ p , w = P^ZlPo, 

Il viendra : 

^=F[p 2 ( Pl -p )-p(pf-p§) + Pi=iPi]- a ^+2^p- r . 

Pour que B soit minimum, il faut d'après les règles du calcul des variations que 
hB s'annule toutes les fois que J est nul et pour cela il faut que l'on ait : 

d % V 

K étant une constante qu'il reste à déterminer. 
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Voyons comment le problème pourra être traité. 
Posons d'abord 

nos équations deviendront : 

,5"= V[2a\ + -^-] — a'a* + 2(3'a — y', 
=f Vda , P'=f Vada , y'=f Vo*d<j , 



ou: , 

a' 



les intégrales étant prises entre les limites — o et. + a . 
Posons maintenant a = a t, 

Vdt, (3"=J Vtdt, y"=J Vfdt. (13) 

L'équation devient: 

L'équation (14) contient encore quatre indéterminées, a", (3", y" et K 1 . Si nous 
l'intégrons nous trouverons : 

V=a!'V 1 + $Ï'Y % + y"V 3 + 5"F 4 + e"F B , 

T^lj T^a» 1^3» T 7 ^ et T^ étant des fonctions entièrement connues de t et de K' pen- 
dant que h" et s" sont deux constantes d'intégration. 

Pour achever de connaître Fil nous restera à déterminer les six constantes 
a", (3", y", è", s" et K '. Pour cela nous avons six équations ; à savoir, les trois 
équations (13), l'équation A = 1 et les deux relations 

— = pour p = p , —■= pour p = p x . 

Ces six équations ne suffisent pas toutefois pour déterminer complètement ces 
six constantes et en particulier EJ. On trouve pour K 1 une infinité de valeurs 
positives. Nous prendrons la plus petite de ces valeurs que j'appellerai K . 

Je n'ai pas besoin pour mon objet de calculer effectivement K ; il me suffit 
de faire observer que c'est une constante numérique. 

Il vient alors : ^_ ^ _ 4iT„ 



4 (pi—pvf 
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Il nous reste à chercher le minimum de -j- correspondant à cette valeur de K. 

XL 

Nous trouvons : 
ou B = KfdoV[vf{ 9 - p'fdp' -/F(p - o'fdp'] 

=^ffy(v-v')( P - P 7d 9 d^. 

Or l'intégrale du second membre ne doit pas changer quand on permute p et p' ; 
on a donc aussi 

B = - Kff F ( F-T')(p — p') 2 d 9 d 9 \ 
d'où B = Kjj { v __ m _ p? ^ = ^ ^ 



Ainsi le minimum de — r- est égal à 



2f SLK, 



2 (/>!-/>o) 5 

Pour une fonction V quelconque on aura donc : 



a >i^rr < 16) 



Soit A, la plus grande distance de deux points de la surface du corps solide envi- 
sagé on aura b 2K 



Il est à remarquer que si je n'avais pas voulu indiquer sommairement la manière 
de calculer la constante numérique K Q , j'aurais pu arriver à la formule (15) en 
quelques lignes par de simples considérations d'homogénéité. 

Il suit de là que l'expression (11) est toujours plus grande que 

6K W 

Par conséquent pour un solide convexe quelconque on a : 

6K a W 



h> 



iùi 



5 i 



K Q désignant une constante numérique, W le volume du corps, et % la plus grande 
distance de deux points de la surface du corps. 
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Cela posé passons à un solide quelconque ; on peut le décomposer en n — 1 
solides partiels convexes. On calculera pour chacun de ces solides le rapport 

W • W 

—g- ; imaginons que pour tous ces solides — r soit plus grand que a ; on aura pour 

A A 

le solide total 6iT n 

Or nous pouvons prendre n assez grand et choisir nos n — 1 solides partiels 
de telle sorte que la quantité que nous venons d'appeler a soit aussi grande que 
l'on veut. 

Donc h n sera également aussi grand que l'on voudra. 

Donc pour un solide quelconque k n croît indéfiniment avec n. 

Nous avons démontré ce théorème dans le cas. où #, = 0; cel'a doit suffire; 
car k n est croissant avec h ; le théorème peut donc être regardé comme démontré 
pour toutes les Valeurs positives de h. 

Je ne veux pas quitter ce sujet sans avoir indiqué un moyen de calculer une 
limite supérieure de h n . 

Posons F = ouFl + a % F % + + a n F n , 

F 1; F%, . . . . , F n étant des fonctions données et a x , a 2 , . . . . , a n des coefficients 
indéterminés. 

Les deux intégrales : 

*=*/**+/[£)•+ (£)'+(£)>. *=/"> 

seront des formes quadratiques dépendant des n paramètres a x , a 2 , . . . . , a n . 
Formons la forme quadratique : 

où Â est un nouveau coefficient indéterminé. 

Ecrivons que le discriminant de la forme B — %A est nul ; nous obtiendrons 
une équation algébrique d'ordre n en a ; il est aisé d'établir que cette équation 
a toutes ces racines réelles (parce que les deux formes B et A sont définies posi- 
tives ; il suffirait d'ailleurs que l'une d'elles le fût. Mais de ce que les deux 
formes sont toutes deux définies positives, il résulte que les racines sont non- 
seulement réelles, mais positives). 

Soient Âj, Â 2 , , Â„ 

34 
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ces n racines rangées par ordre de grandeur croissante ; je dis qu'on aura : 

^1 !> \ y ^2 ^> ^2 ) ^3 ^> #3 ! • • • • > %>n ^> *>n • 

En effet la plus petite valeur que puisse prendre le rapport -j- quand on fait 
varier les a doit être plus grande que \ . Or cette plus petite valeur est \ ; on 
a donc X 1 >ft l . 

En vertu de la théorie des formes quadratiques, la forme A peut être décomposée 

en n carrés et s'écrire : 

A=ft + ft + ....+pi, 

fi lt /8 g j3 B étant des fonctions linéaires de a lf a 2 , a„. 

De plus cette décomposition peut être faite de telle sorte que 

Introduisons maintenant les conditions : 

fFUidrt =fFUzck = =fFU v dv = . (16) 

Le minimum de -j-, en tenant compte des conditions (16), devra être plus grand 

que h p + 1 . 

Cherchons ce minimum. Les conditions (16) sont linéaires par rapport aux 
(3 de sorte qu'on peut les écrire : 

l«n/?i + ^i2& + ••••+ /W?» = 0, 

^2lPl + ^22^2 + ....+ [hnPn = 0, /jg ^j s 

Soit à le minimum cherché du rapport —r- . On voit par des calculs qu'il est 

inutile de reproduire ici que ce minimum qui nous est donné par l'équation 
suivante, dont je vais expliquer la signification : 

XH*Il(X — ù i ) = Q. (17) 

EL (X — ht) sera le produit de n — p binômes tels que % — ■. \. Supprimons dans 
les équations (16 bis) tous les (3 t qui ont même indice que les X { qui entrent dans 
le produit II (/l — y^) ; nous appellerons H le déterminant des équations linéaires 
ainsi obtenues. Enfin la sommation indiquée par le signe 2 est étendue à toutes 
les combinaisons des n quantités \ prises n — p à n — p. 
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On voit qu'en substituant dans l'équation (17) successivement 

— oo, a, 1( Xi, . . . . , K, + °° 

on obtient au moins p changements de signe. Par conséquent la plus petite 
racine de l'équation (17) sera au plus égale à \ +1 . On a donc : 

%>p + l>k p + 1 . C.Q.F.D. 

Je termine ce paragraphe en donnant une nouvelle manière de calculer une 
limite supérieure de h % pour h == . 

Soient/, g, h trois fonctions quelconques assujetties à la condition suivante : 
A la surface du corps on aura : 

af+fig + yh=0, (18) 

a , fi , y désignant les trois cosinus directeurs de la normale à cette surface. 
De plus nous assujettirons nos trois fonctions à la condition 

A=f(f + Î ?+Jf)dv=l. (19) 

I»,e *=/[(£)■+(£. )'+(!)> 

aura évidemment un minimum. Cherchons ce minimum. Posons : 

g_ d f ■ % ■ dh 
dx dy dz 

Nous trouverons, par le calcul des variations : 

aèf+ fièg + yèh = (à la surface du corps) (18 bis) 

et \ U =f(fèf + gèg + hèh) ck=0, 

2 *f \ dx dy dz J 

L'intégration par parties donne : 

^8B=f6(aèf+fièg + yèh)â«-f(§èf+^Sg+^8h)d* t 
ou en vertu (18 bis) 
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Pour que SB soit nul toutes les fois que è A est nul, il faut donc que l'on ait : 

7 /• . dO , , dd „ 7 T , dd , „ „ 

A; étant une constante qu'il reste à déterminer. 

Des équations (20) on tire par différentiation et addition : 

A0 + k$ = 

et l'équation (18) devient : d0_ n 

dn 

Cela montre que est l'une des fonctions TJ que nous avons définies plus haut ; 
ce ne peut être que Z7 2 ; on a h = h % , et pour $ = TJ % on trouve : 

B f&ck JZfufdrt 

Pour des fonctions/, <?, A quelconques on aura donc 

d'où la règle suivante. 

On prend trois fonctions quelconques f, g, h assujetties à la condition (18). 

La condition (19) devient inutile dès qu'on considère le rapport — r- J . Le 

rapport /Vj/ dg 6Jk N* 

t/ vas ay efe / 

f(f + f + h*)dx 
est plus grand que &j . 

§3. — iois cZt* Refroidissement. 

Soit F la température d'un point du corps solide ; V sera une fonction de 
x, y, z et de t. Cette température devra satisfaire aux deux équations suivantes : 

dV 
A l'intérieur du corps: (1) — — =a 2 AF. 

dV 
A la surface du corps : (2) — — |- h F= . 
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La température est donnée arbitrairement à l'époque initiale t = . Il peut 
donc se faire qu'à cette époque t = 0, l'équation (2) ne soit pas satisfaite ; mais 
elle devra l'être pour toute époque postérieure t > . C'est là une première 
anomalie qui vaut la peine d'être remarquée. 

En voici une seconde : Y ne peut pas en général être développée suivant 
les puissances croissantes de t. Supposons en effet que cela soit possible; 
qu'arrivera-t-il ? Soit Y la valeur de Fpour jf = . On aura pour t = 0: 



dV „. „ cPV 2A dV 



= « 2 AF 0l -^ = a 2 A — = a 4 AA7 



ou 



cPV 

d? 



= «'A 8 7 



en convenant de poser : 
On aura ainsi en général : 



A n V=A(A n ~ 1 Y). 
_L — n 2n A n V • 



dl n 

de sorte que si le développement était possible, il viendrait : 

V= V + aHAY + ^A*V + J^A 3 F + . . . . 

Il résulterait de là que la température en un point donné et à un instant donné 
ne dépendrait plus que de la valeur de Y et de ses dérivées en ce point. La 
forme du corps solide n'interviendrait en aucune façon. Cela est absurde. 

Pour mieux faire comprendre ces anomalies, nous allons envisager un cas 
particulier. Imaginons que le solide devienne un mur indéfini compris entre les 
deux plans x= ±7t. 

Supposons que les deux plans x = db it qui limitent le mur soient impermé- 
ables à la chaleur ce qui revient à supposer h — . 

Supposons de plus qu'à l'époque t = ; la température initiale V ne dépende 
que de x et ne change pas quand on change x en — x . 

Ces propriétés subsisteront évidemment à une époque quelconque. Y sera 
fonction de x et de t seulement et ne changera pas quand on changera x en — x. 

Dans ces conditions on peut poser : 

V= 2<p m (t) cos mx . 
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Nous aurons alors A Tr <&V „ « ,. 

-T- = 2<f4(/)cosma;. 
Si dans l'équation dF a „ 

nous faisons a 2 = 1 pour simplifier il vient : 

4>U*) + «ty.(*) = 
d'où 4, m («) = ^- ma< 

et enfin F= 2 J. m e~ ro2|! cos mec . 

Donnons-nous la température F à l'instant initial t = . F pourra toujours 
être développée par la série de Fourier sous la forme : 

F = 2^. m cos mx . 

La série du 2 a membre ~ZA m cos mx est toujours convergente, mais la conver- 
gence peut n'être pas absolue. 

En vertu d'un théorème d'Abel, si l'on prend 

V— 2A m e- mn cos mx (3) 

on aura quand t tendra vers 

lim V= V Q . 

L'équation (3) nous fournit donc la solution du problème. 

Il semble que la condition A=0 n'ait joué aucun rôle dans ce calcul; ce 
n'est là qu'une apparence à la quelle il ne faut pas se tromper. 

Nous pouvons, il est vrai, dans tous les cas possibles, développer Y par la 

série de Fourier, et écrire : 

V= 2$ m (t) cos mec. 

Mais pour que nous ayons le droit d'en conclure 

j-j- = — 2wi 8 $ m (t) cos mx 

il faut que la série (et d'ailleurs il suffit) 

2»n 2 4v (t) cos mx 
soit convergente. 
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Pour cela il suffit que l'on puisse trouver un nombre K tel que 

l"»**«(OI <-^- 

l*.WI<-§. (4) 

D'après un théorème que j'ai démontré dans le Tome III du Bulletin Astro- 
nomique (Sur un moyen d'augmenter la convergence des séries trigonométriques) 
la condition (4) équivaut à la suivante ; la fonction représentée par la série 

/(ce) = 2<p m cos rnx 

devra être continue ainsi que ses trois premières dérivées. Or cette fonction est 
égale à F entre les limites — n et + n; si on est en dehors de ces limites on a : 

f(x)= V(x+2jp7t), 

p étant un entier positif ou négatif choisi de telle sorte que ce + 2p7t soit compris 
entre — n et + n . 

Comme V est continue ainsi que toutes ses dérivées, il ne peut y avoir de 
discontinuité qu'aux deux limites x = ± n . Si donc nous désignons pour un 
instant par V, V", etc. les dérivées successives de V par rapport à ce, on devra 
avoir: y (n) - V (-n), (5) 

V (7t)= V (-7t), (6) 

V" (n)= Y" (— n), (7) 

V'"{n)=V"i(—7t). (8) 

La fonction Fêtant paire les conditions (5) et (7) sont remplies d'elles-mêmes. 
D'autre part V et V" sont des fonctions impaires de sorte qu'on doit avoir : 

V (7t) = — V (- n) , V»> (n) = - V" (— 7t) 
et que les conditions (6) et (8) peuvent s'écrire : 

V'(7t)= F"(7l)=0, 

c'est-à-dire que pour x = n on devra avoir : 

dv__ârv_ ,, 

dx dx* ' ^ ' 



Si h = , on doit avoir pour x = n : 

dV _ dV 
dn dx 



= 0. 
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Cette condition ayant lieu, quel que soit t on aura : 

d?V_d?V_ 
dtdx dx 3 

Les conditions (9) sont donc remplies, et notre calcul est légitime, mais seulement 
dans le cas de h = . 

Cela posé considérons la série 

V = %A m cos mx. 

En général la série : Xm 2 A m cos mx 

ne sera pas convergente de sorte qu'on ne pourra pas écrire : 

cPVo 
dx* 



AF = -j-^ = — 2rn 2 A m cos mx . 



Il en résulte qu'on n'aura pas en général : 

lim AV= AV quand t tend vers 

et qu'on n'aura pas non plus : 

lim -— = AV quand t tend vers 0. 

C'est ce qui explique pourquoi le développement suivant les puissances de t 
est généralement impossible. 

Revenons maintenant au cas général : 

On a pour t > dV , „ 

dn +^- u > 



ou en différentiant par rapport à t: 



d?V , _ dV n 



dtdn dt 

ou: dAV 



dn 



+ hAV=0. 



En différentiant p fois on trouverait de même 

dA p V 



dn 



+ hA p V=0. 



Pour que le développement suivant les puissances de t soit possible, il faut 
évidemment que l'on ait : 

^ + hV =0, d ^ + hA*>V (j =0. (p = l, 2,....,adinf.) 
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Ces conditions, en nombre infini, sont nécessaires ; j'ignore si elles sont suffisantes 
quoique cela puisse sembler probable. Je n'insiste d'ailleurs sur tous ces points 
que pour mieux montrer avec quelles précautions il faut toucher aux équations 
aux dérivées partielles. 

Passons maintenant à l'exposé des lois générales du refroidissement. 

Considérons d'abord l'intégrale suivante : 

A=fV*àW, 
Je dis que cette intégrale ira toujours en diminuant ; nous trouvons en effet : 

Il vient ensuite : 

ou en vertu de l'équation (2) 

/VAVck=—hf V 2 da —f$ (j^dx =-B<0 

et par conséquent : ^ = - 2alB < . 0. Q. F. n. 

ai 

Je dis maintenant que si h n'est pas nul, / V*dr tend vers quand t croît 

indéfiniment. On a en effet B_ 

A 

d'où dA 



>h 



A di <- 2 ^ 
ou en appelant A la valeur de A pour t = t , 

A<A & e- Zak > t , 

Si h n'est pas nul, \ n'est pas nul non plus, et l'on a : 

lim A = pour t = oo c. Q. F. d. 

Je dis maintenant que le rapport — va constamment en diminuant : 
Il vient en effet : d / B \ AdB—BdA 



dt \A ) 



AHt 
35 
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de sorte qu'il s'agit de démontrer l'inégalité suivante : 

. dB „ d A . ,„ 

A — B—r <0. 

dt dt 

Nous devons donc d'abord calculer ~-=~- , il viendra : 



dB 



— tf™=-f%*y*-fY!Ç*. 



dt 

L'équation (1) nous permet donc d'écrire 

ldB 



aHt 



= -f(AVfdx—fvA* 



Vdt. 



Or si V et JJ sont deux fonctions satisfaisant à la surface du corps à l'équation 
(2) le théorème de Green nous apprend que 

f(VAU— UAV)dt-0. 

Mais F et A F satisfont à l'équation (2). On a donc 

f VA 2 Vdt =/( A Vfdt 



et par conséquent : dB^ 



/(AT) 2 



= — 2a? / (AV) 2 dt 



dt 

ce qui montre déjà que B est décroissant. 

Nous pouvons écrire (en appelant aW un élément de volume du corps autre 
que dt et désignant par V la valeur de V au centre de l'élément du 1 ) 

^ = 2a 2 / VA Vdt' ; ^§--~ 2a 2 /(A V'fdt 1 
et par conséquent : 

A ~ - B ~ = - 2a 2 [/ Vdtf(A V'fdt 1 -f VA Vdtf VA Vdt'"] 
ou A~~B~z=- 2a*ff[( VA Vf - VVA VA V] dtdt*. 

Comme rien ne distingue dt de dt', nous pouvons écrire également: 

A ^dT~ B ^dT-~ 2a *fftt V ' A F) * ~~ FFA FA V 1 dtch '' 
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d'où ajoutant les deux valeurs de A -j B —5— et divisant par 2, 



dt dt 



*d*d*' 



d'où . dB „ dA . n 

A~ J5-— <0. C.Q.F.D. 

dt dt 

Nous avons vu plus haut que l'on a : 

A < ArfT**** 

cette inégalité peut dans certains cas être remplacée par une autre. Supposons 

que l'on ait: /» /» /» „ 

J VU x âx —j VU % dr = —J VU n _ x àW = 

nous aurons d'après définition même des quantités h n et U n : 

il viendra donc : dA = _ ^ ^ 

dt 
et A<A e~ Za ^K 

Etudions maintenant les variations de l'intégrale : 

J n = J*VU n ck. 

Ilvient: ^=a>fu n AVd*. 

Le théorème de G-reen donne : 

En vertu des égalités 

le premier membre est nul ; on a donc : 

/ U n A Vck =f VA U n d* = - h n f VU n dt = - hj n 
d'où dJ^_^ ,, r 

et J n = Jle-^\ 

Jl représentant la valeur de l'intégrale J n pour £ = . 
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Etudions encore les variations de l'intégrale 

H=fv i V i dt 

oxx V 1 représente la température à l'instant t et V % la température à l'instant 
h — t. Il vient: 

^ = a 3 /(F 8 AF 1 -F 1 AF,) ( fo 
or en vertu du théorème de Green et des équations 

le second membre est nul. Donc H est une constante qui ne dépend que de h . 
Si nous faisons . -, . h 

il vient : V 1 = V z 

et H=fvi&t>0. 

Si donc V x et Y % représentent les températures à deux instants quelconques, 
l'intégrale /» „ „ 

sera toujours positive. 

Nous avons vu plus haut que l'intégrale 

considérée comme fonction du temps va toujours en décroissant. 

A 

2 



Donc H qui si l'on fait t = — - se réduit à 



fvfck 
sera une fonction décroissante de h. 

Or nous trouvons : dH ./»„ Al77 y „ 

Nous trouvons de même 

f V % A V x èe =f V,A F 2 <fe < o . 
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Si donc V x et V % représentent les températures à deux instants quelconques, 
l'intégrale /» „ 

sera toujours négative. 

Arrivons maintenant au problème principal ; étant donnée la température 
au temps t = , trouver la température à un instant quelconque. 

Soit V la température à l'instant t = . 

La solution classique consisterait en ceci : 

Développer V en une série de la forme suivante : 

V = A 1 U 1 + A i U,+ ....+A n U n + ...., 

les A étant des constantes. 

On aura ensuite à un instant quelconque : 

V= A x e~ a% ^U x + Ajr^Uz + + A n e~ a ' M U n + 

Cette solution est subordonnée à la possibilité du développement, et c'est cette 
possibilité que nous ne sommes pas encore en mesure de démontrer d'une 
manière générale. 

Voici toutefois ce que nous pouvons dire. 

Soient A 1} A 2 , A n des coefficients quelconques; posons: 

V = A 1 U 1 + A,U z +....+A n U n + B 

et proposons-nous de déterminer les coefficients A de telle fapon que l'erreur 
moyenne commise soit minimum. 

Nous prendrons, à l'exemple de M. Tchebicheff, pour mesure de l'erreur 
moyenne commise l'intégrale suivante : 

S =fl%d*. 
Cherchons donc le minimum de l'intégrale 

f(V -A 1 U 1 -A z U z -....-A n U n fdr. 

Cette intégrale sera minimum quand on aura : 

fu p (V -A p U P )d* = O 
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ou (puisque nous avons par définition 

quand on aura : A p = J£. 

Nous sommes donc conduits à écrire : 

Vo^JÏU. + JIU, + ....+ J° n U n + B . 

Il résulte de là que l'erreur moyenne commise S va en diminuant quand n 
augmente, mais non que cette erreur moyenne tende vers quand n croît au 
delà de toute limite. D'ailleurs S pourrait tendre vers sans que B tendît 
vers 0. 

Remplaçons toutefois V Q par sa valeur approchée 

p = n 

p=i 

Nous en déduirons £=^ ^-^ 

p=i 
nous rappelons que nous avons posé 

f VTJ p ck = J p = jy***. 

Posons donc : V= J-JJ X +J % U % +.... + J n U n + B, 

et prenons pour mesure de l'erreur moyenne commise : 

S=fB*ck, 

je me propose de démontrer que l'on peut prendre n assez grand pour que 
l'erreur moyenne S commise sur la température à un instant donné soit aussi 
petite qu'on le veut. 

En effet B satisfait aux équations 

~ = a*AB, ^ + hB=0. 
at an 

Si donc la température initiale était B , la température à un instant ultérieur 
serait représentée par B. 
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De plus on a, comme il est aisé de le vérifier: 

fs Utân =f% V z ck — =/b UJ* = . 

Donc d'après un lemme que nous avons démontré plus haut, on aura : 

S<S e- a, ^ { . 

Or quand n croît au delà de toute limite, S décroît (sans que nous sachions 
s'il tend vers 0), k n+1 croît au delà de toute limite et l'exponentielle e - o2J: »+i* tend 
vers 0. Donc S tend vers 0. o. q. f. d. 

Nous avons donc démontré que l'erreur moyenne S tend vers , mais non 
que B tend vers 0. Cela peut toutefois nous suffire pour le moment. En effet 
comment pourrait-il arriver que S tendît vers sans qu'il en fût de même de Bt 
Il faudrait pour cela que la valeur de B subît des oscillations d'autant plus rapides 
que n serait plus grand, de telle fapon que pour n très grand, B prendrait en des 
points très rapprochés des valeurs très différentes. Aucun physicien ne doutera 
que si un pareil état de choses existait à l'instant initial, il ne saurait subsister. 
C'est ce qui m'engage à me contenter pour le moment des considérations qui 
précèdent. 

Je terminerai ce paragraphe par la remarque suivante : 

Si V est partout positif, V sera aussi positif pour toutes les valeurs de t et 
pour tous les points du corps. Or quand t croît indéfiniment, le rapport : 

V 



tend vers l'unité. 

Donc Z7i doit être une fonction qui est positive en tous les points du corps. 

L'égalité fu 1 U n dx = (»>1) 

montre que la fonction TJ X est la seule des fonctions U n qui jouisse de cette pro- 
priété. 

§4. — Propriétés des Fonctions U n . 
Reprenons la fonction U n définie par les équations : 

AU n + k n U n =0, ^ + hU n = Q, fïïlék=l. 
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Ou en supprimant les indices qui nous sont inutiles pour le moment : 

Soit T une fonction satisfaisant comme U à l'équation : 

AT+kT=Q. 



Si T est finie et continue ainsi que ses dérivées à l'intérieur du corps, on aura 





U^-T^d»: 



dn dn / 



/( 



ou 

*' "\dn """'"V 

les intégrales étant étendues à la surface du corps. 

Supposons maintenant que la fonction T ne soit plus finie à l'intérieur du 
corps qu'elle devienne infinie au point (x , y , z ) situé à l'intérieur du corps; 

mais de telle façon que la différence T (où r désigne la distance des deux 

points x, y, z et x , y , z ) reste finie ainsi que ses dérivées. 

On aura alors : f(u^- T^-) do = - 4*17», 

v \ dn dn / 

U° désignant la valeur de U au point x , y , z . C'est ce que j'ai déjà exposé à 
propos de la Diffraction dans ma Théorie Mathématique de la Lumière. 

In vient donc : fu(^ + AT) do = - AnU\ 
Soit maintenant : a = s/ h 

et m ^ 

1 ~~~T' 
r désignant encore la distance du point mobile x, y, z au point fixe x , y , z . 



On aura alors 



fu(^- + AT\do=0 ou — 4nU°, 



selon que le point x , y , z Q est extérieur ou intérieur au corps. 
On a d'ailleurs dans ce cas : 



= cos4— (^-—J, 



dT .... , é 
dn 
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4- étant l'angle de la normale à la surface du corps au point x, y, z avec la droite 
qui joint ce point au point fixe x , y , z . 

Nous poserons pour abréger 

dT 
H=~+hT 

dn 

et nous regarderons H soit comme fonction de r et de 

r cos 4 1 == £ 

soit comme fonction de x, y, z et de x , y , z . 
On aura alors : 

r cos ^ = à (ar — x) + fi (y — y) + v(z — z), 

h, [i , v désignant les trois cosinus directeurs de la normale à la surface du corps. 
On en déduit : dff _ dH x — x dH 

dx dr r dç 

On a ensuite si le point x , y , z est intérieur au corps : 

471 v ' dx 47t «/ dx 

Cela va nous permettre de trouver une limite supérieure de U et de ses dérivées. 
Soit en effet A la plus grande valeur que puisse prendre \U\ à la surface 
du corps, on aura évidemment : 

JTT A /» I J m 

do. 



i^k-^/Vk», 


dU 
dx 


< 


A P\dH 
An J 1 dx 


Les deux intégrales qui entrent dans ces deux inégalités 


f\H\doet t 


f\i 


IH 

lx 


do 



se calculent aisément quand on connaît la forme de la surface du corps et le 
nombre positif 7e. Elles ne dépendent que de x , y , z . 

Quant au coefficient A, nous n'avons jusqu'à présent aucun moyen de le 
déterminer. 

Commençons par étudier les variations de la première intégrale / |j5~|e?a. 

Cette intégrale est évidemment finie tant que le point x , y , z reste inté- 
rieur au corps. On a en effet : 

36 l ff l<7 + 7 + 7 
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de sorte qu'il viendra : 

S désignant la surface totale du corps et p la plus courte distance du point 
œ o> Vo > a o à cette surface. Je dis maintenant que notre intégrale tendra encore 
vers une limite finie quand le point x , y a , z se rapprochera indéfiniment de cette 
surface. 



Posons en effet : 






H x sera une fonction qui ne deviendra pas infinie même quand r s'annulera. 
Nous trouvons en effet : 

COSi^ 



J2i = 



[taré* — e Ur + 1] -\ {é ar — 1) , 



e iar_ 

r 


• 1 


2 


ar 
sm — 


-<a; 




r 



iare iar — e iar + 1 



<4a 8 



et enfin 
d'où : 



\H 1 \<ha + 4a z , 

Il est aisé de voir que quand le point x , y , z se rapproche indéfiniment de la 
surface du corps, les intégrales / L_ — zA e t l'intégrale / — tendent vers des 

limites finies. 

Si d'abord le corps est convexe, de telle façon qu'une droite ne puisse couper 
sa surface en plus de deux points ; cos 4 1 est positif et l'on a : 

P Gos ^dco==f^dc> = 4n. 

Si le corps n'est pas convexe, et qu'une droite puisse rencontrer sa surface en n 

points on aura: /*|cos^| , . . , ,. 

J -1 — -^-^i- cto < 4 (n — 1) 7ï 

car l'intégrale s'obtient en décrivant du point x , y , z comme centre une sphère 
de rayon 1 et en faisant la perspective de la surface du corps sur la surface de 
cette sphère, le centre de la dite sphère étant pris comme centre de la perspective. 
Un point de la sphère sera alors au plus la perspective de w — 1 points de la 



Partielles de la Physique Mathématique. 277 

surface du corps; de sorte que la somme arithmétique des aires des perspectives 
des divers éléments de cette surface sera au plus égale à n — 1 fois la surface de 

la sphère. Or cette somme arithmétique est précisément l'intégrale, / -! — ^J da. 
La somme algébrique serait l'intégrale f °\ da. 

Ainsi / -1 — -g-^J- da tend vers une limite finie ; il reste à démontrer qu'il en 

est de même de f* da 

J r 

Or cette intégrale représente le potentiel d'une couche, uniforme de matière 
répandue à la surface du corps, et l'on sait que ce potentiel est fini. 
Si nous passons maintenant à l'inégalité : 



dUo <A f\dH 
dxo J I dx 



da, 



elle nous permet de démontrer qu'à l'intérieur du corps les dérivées premières (et 
on le démontrerait de la même fapon pour les dérivées d'ordre supérieure) de la 
fonction U restent finies; mais elle ne nous permet pas de voir si ces dérivées 
tendent vers une limite finie quand le point x q , y , z se rapproche indéfiniment 
de la surface du corps. 

Nous allons maintenant chercher à obtenir une limite supérieure du coeffi- 
cient A . 

Pour cela il nous faut démontrer que U est une fonction continue. 

Cela est évident pour les points situés à l'intérieur du corps puisque nous 
venons de voir qu'en ces points on peut trouver une limite supérieure des dérivées 
de U. 

Il reste à démontrer que U est encore une fonction continue sur la surface 
du corps, et pour cela il nous faut une expression de U quand le point x , y , z 
est sur cette surface. 

Nous avons trouvé plus haut 

fjJHda = ou — 47tZ7 
selon que le point x , y , z est extérieur ou intérieur du corps. 
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En vertu des mêmes principes, on aura 

fuHd(o = — 2nU 

si le point x , y , z est sur la surface même du corps. 

Si x' , y(,, z' est un autre point situé également sur la surface du corps et si 
H' et Z7o sont deux fonctions formées avec le point x(,, y!,, Zq comme H et U \e 
sont avec le point x , y , z ; il viendra : 

f U{H— H') da = 27t(Ul ) — U ) 

d ' où: \Ui—U \<Af\ff—H'\d a . 
Il nous reste à établir que : Ç , jt jt, i j 

tend vers quand le point x' , y' Q , z' 6 se rapproche indéfiniment du point x , y , z . 
Remarquons d'abord que nous pouvons poser : 

a- = ^._ ««* + «, 

.Hi a même signification que plus haut; H{ est formé avec le point x' , y' , z' 
comme jE^ avec le point x , y , z ; r* désigne la distance du point x' , yô, z' Q au 
point x, y, z et l'angle 4 1 ' est défini avec le point x' , yl, z' Q comme l'angle 4> avec 
le point x , y , z . 
On a donc : 

/|j-ir|<b<»/|i— ^|&+/|~4_!^4!|A+/|ii;-flïl&. 

Il suffit donc de démontrer que les trois intégrales 

/lr-7-1*- /|°-^-^|*. /l«-iS|4. 

tendent vers quand les deux points se rapprochent indéfiniment. Cela est 
évident pour la troisième; car H x est une fonction continue et finie de x , y et z . 

Démontrons-le maintenant pour la première. 

Décomposons la surface S du corps en deux régions a et <r' ; et supposons 
que les deux points x , y , z et xb, yl, z& appartiennent à la région cr. 
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Je dis que je pourrai prendre ces deux points assez rapprochés pour que 



s I r 



da'- 



dot + / , -r 



da<C s. 



Tout d'abord nous savons que les deux intégrales 



Pda pda 



sont finies et déterminées. Il en résulte que nous pourrons prendre la région <j 

assez petite pour que : 

Jf*doi . s Pda . s 

quelle que soit la position du point x^, y^, z' dans cette région a d'où : 

2e 



f.\ 






C?6>< 



La région a doit désormais être regardée comme entièrement déterminée, 
mais nous pouvons encore faire varier dans cette région le point x'q, yi>, z' . 

Si maintenant r 1 désigne la distance du point Xq, yh, z[ à un point ce, y, z 

quelconque de la région a' ; —j- sera une fonction finie et continue de xl, y' , Zq et 

cette fonction tendra uniformément vers — quand Xq, y{, z[ tendront vers x , y , z . 

Cela sera vrai tant que le point x, y, z restera sur la région a', puisque les points 
x o> i/01 z o e * x o> yô> z o n'appartiennent pas à cette région. 

On pourra donc prendre le point Xq, y' , z' assez voisin du point x , y , z 
pour que 



V «r 7 I f y* 



et par conséquent : 



r 



do< — 



da<C s- 



0. Q. F. D. 



On établirait de la même manière que 



/' I cos 
im J \-? 



lim 



4> cos 4' 



On a donc : 



c&j<0. 

lim|£7- -£7- '| = O 



limf\H—H'\da = 
ce qui signifie que la fonction U est continue à la surface du corps. 
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Cette démonstration ne montre toutefois qu'une chose, c'est que TJl tend 
vers TJ , quand le point x , y a , z appartient à la surface du corps et que le point 
£Co, 2/0» 3 o s ' en rapproche en restant lui-même sur la surface du corps. Elle ne 
nous apprend rien sur ce qui arrive quand le point x' , yl, z' a est intérieur au 
corps et se rapproche de x , y , z soit normalement, soit obliquement à la surface 
du corps. On pourrait observer toutefois que 

an 

étant fini, TJl doit tendre vers TJ quand la droite qui joint les deux points est 
normale à la surface et il serait aisé d'en conclure, par un petit raisonnement 
très simple, qu'il en est encore de même quand cette droite est oblique. 

Mais cela ne saurait nous suffire, parce que nous avons besoin pour notre 
objet d'assigner une limite supérieure de | TJl — TJ e \. 

Si comme nous le supposons le point x , y , z^est sur la surface du corps et 
le point x' 0l y' , z' à l'intérieur, on aura : 

J*HTJda — 2nTJ =: — 4* TJ , 

fH'TJdo — — AnUl 

d ' où f(H - H') TJdu —2nTJ Q =—4n(TJ — VI) , 

Nous avons donc : 

- 4n ( TJ - TJl) = hf(±- — £.) OS» 

+f(ff 1 - HQ Udo +f(^L - ^) Uda> - 2nU . 

Il nous faut démontrer que les trois intégrales du 2 d membre tendent respective- 
ment vers 0, et 27iTJ Q quand le point x' , y' , z^ se rapproche indéfiniment du 
point a- , y , z . 

Or on a: \f(H x - H() Udu < Af\H l -H{\ 

et on verrait comme plus haut que 

limy \H l — 3{\â0zsQ. 



\do 
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Il n'y aurait rien à changer à la démonstration que nous avons donnée dans le 
cas qui précède. 

De même on a : 



et 



lt/\y r' / ** I r ¥ 

/Il 11 
-r du=0. 



d(ù 



Ici encore il n'y a rien à changer à la démonstration que nous avons donnée dans 
le cas précédent. 

Il reste à établir que : 



Bm/C^-SJ*)!»-**. 



Décomposons la surface S du corps en deux régions a et a' et supposons encore 
que le point œ , y , z soit situé dans la région a; il viendra : 

fft? - *£) **• =£&? - ^) ud ° 

Je dis que je puis prendre le point x(,, y' , z' assez voisin de x , y Q , %pour que la 
différence de 2nU () et du premier membre de l'identité (2) soit plus petite en 
valeur absolue que s. 

En premier lieu, la fonction U étant continue à la surface du corps, nous 
pourrons prendre la région <r assez petite pour que sur cette région on ait : 

\u-u \<ç, 

Ç étant une quantité aussi petite que l'on veut. Il vient alors : 



|/(^-^)(^.H<î/|^-^ 



do. 



Or nous venons de trouver 



/| cos 
r 



COS'vf' 



da < 4 (n — 1) n, 



n étant le nombre maximum des intersections d'une droite avec la surface du 
corps. 
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Il vient donc : 

\f.Q^ -^){u-u)^\< H »-^<^. 

Car, Ç étant aussi petit que je veux, je puis toujours prendre 

s ^24(«— i)n' 

La région <r doit être regardée maintenant comme entièrement déterminée, mais 
je puis encore faire varier le point «é, y(,, se- 

Nous voyons d'abord que le point x , y , z ne faisant pas partie de la région 

c'. on a: /*|cos^ cos^ , 

de sorte qu'on peut prendre le point xô, yô, zô assez voisin du point x Q , y Q , z pour 



que: 



cos -4/ cos 4> 



do < 



SA 



et par conséquent pour que : 

I /VcoH/ _ coMA & I e 
I t/<r'V r 1 r* J 3 



D'autre part l'intégrale 



' cos j> ; 



C?6) 



réprésente l'angle solide sous lequel on voit le contour de la région <r du point 
»o> yô, 2o» et l'intégrale 

réprésente l'angle solide sous lequel on voit le même contour du point x , y , z . 
Il en résulte que : 

lim J^ 3 da = / -~ da + 2ti. 

On peut donc prendre le point x^, y' , z[ assez voisin de x Q ,y , z pour que 
I /Vcos i|/ cos 4-\ ^ I • e 
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Il vient alors I £(^_ _ W) Uda _ 2yt Uo 



^ 3 



<C e. C. Q. F. D. 

Voici le résumé de cette longue discussion. 

Soient x , y , 2 ; x^, yl, Zq deux points situés soit à l'intérieur du corps, soit 
à sa surface. Soient U et Uq les valeurs de U en ces deux points ; on aura : 

\U,-Ul\<A6, 

A étant la plus grande valeur de | U\ à la surface du corps ; quant à 6 ce sera 
une fonction de x , y , Zq, x' Q , y' , z' que l'on pourra déterminer complètement à 
l'aide des considérations qui précèdent quand on connaîtra la forme du corps. 
Ces considérations nous fournissent en effet une limite supérieure de \U — Ul\ 
puisqu'elles nous montrent comment on doit choisir les deux points x , y , z , 
#o\ Vo i z ô pour que | U — U^\ soit plus petit que s. 

Je n'attirerai l'attention que sur deux des propriétés de la fonction $. Elle 
est essentiellement positive et elle tend uniformément vers quand le point 
*o> yL z ô se rapproche indéfiniment du point x , y , z . 

Regardons d'abord le point x , y , z comme fixe et situé sur la surface du 
corps et faisons varier le point x' , y' , z' . Nous pourrons diviser le volume du 
corps en deux régions que nous appellerons E et R' et que nous définirons 
comme il suit : 

Quand le point x' , y' , z' sera dans la région R on aura : 

0<1. 

Quand ce point sera dans la région R' on aura : 

0>1. 

La région R existe certainement et son volume ne peut être nul, puisque est 
très voisin de quand le point x(,, yô, zl est très voisin de x , y , z . 

Si nous supposons en particulier que le point a- 0) y , z soit celui des points 
de la surface du corps où | U\ atteint sa valeur maximum A; on aura : 

tant que le point xô, y' , z' restera à l'intérieur de la région R. 
37 
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Soit dt 1 l'élément de volume du corps dont le centre de gravité est x' , y(,, %[. 
0naura: fufd*> = l, 

l'intégrale étant étendue au corps tout entier ; et par conséquent 

y>w<i, 

l'intégrale étant étendue seulement à la région R . On en déduit : 

Cette inégalité est vraie pourvu que l'on ait choisi pour le point x , y , z celui 
des points de la surface du corps pour lequel 

\U\=A. 

Malheureusement nous ne savons pas quel est celui des points de la surface pour 
lequel cela a lieu. Mais nous pouvons tourner la difficulté de la fapon suivante. 
L'intégrale /» 

peut être calculée dès que l'on connaît la forme du corps et le point x , y , z . 
C'est donc une fonction de œ , y , z . Cette fonction ne peut jamais s'annuler. 
Elle aura donc un minimum M que l'on pourra déterminer dès qu'on connaîtra 
la forme du corps. Il vient ainsi 

A*M<1 
d'où , ^ 1 



Ainsi nous pouvons déterminer une limite supérieure du coefficient A et par 
conséquent une limite supérieure des dérivées d'ordre quelconque de U en un 
point quelconque de l'intérieur du corps. 

§5. — Retour a V Hypothèse moléculaire. 

Dans les raisonnements qui remplissent les trois paragraphes précédents, il 
y a un point faible que j'ai déjà signalé plus haut. 

Après avoir montré qu'une certaine intégrale ne pouvait pas s'annuler, nous 
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en avons conclu que cette intégrale devait avoir un minimum, et nous avons 
déterminé la fonction U qui correspond à ce minimum par le calcul des varia- 
tions. Or cette application n'eût été légitime que si nous avions démontré 
d'avance la continuité de cette fonction U. C'est d'ailleurs la même objection 
qui empêche de regarder comme rigoureuse la démonstration du principe de 
Dirichlet par Riemann. 

Il est vrai que dans le § précédent, nous avons trouvé une limite supérieure 
de la dérivée de cette fonction U; mais, si l'on voulait s'en servir pour justifier 
l'emploi du calcul des variations, on commettrait une pétition de principe ; tout 
au plus ce résultat peut-il mettre sur la voie dans la recherche d'une démonstra- 
tion satisfaisante. 

Il faudrait donc, pour obtenir une théorie analytiquement rigoureuse, 
employer des procédés analogues à ceux qui permettent d'établir le principe de 
Dirichlet et peut être des procédés plus compliqués encore. 

Je ne l'ai pas fait; mais j'ai pensé qu'il était possible d'obtenir une démon- 
stration rigoureuse au point de vue physique de la fapon suivante. Au lieu de 
considérer l'équation différentielle de Fourier en elle-même, rappelons-nous 
quelle est sa signification physique et comment on l'a obtenue. 

On considère un corps solide formé d'un très grand nombre de molécules. 

Soient M lt M, M n 

ces molécules, n est un très grand nombre. Soient 

Vi, V t , ..... V n 

les températures de ces molécules. 

La molécule M t enverra à la molécule M k une quantité de chaleur égale à 

C ik étant un coefficient indépendant des températures, ne dépendant que de la 
distance des deux molécules ; ce coefficient est très petit dès que cette distance 
devient sensible. 

En outre, cette molécule M t rayonnera au dehors une quantité de chaleur 
égale à: q^ 

Ci étant un coefficient qui n'est sensible que pour les molécules superficielles. 
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Si nous choisissons l'unité de chaleur de façon que chacune de nos molécules, 
que nous supposons toutes pareilles entre elles ait pour chaleur spécifique l'unité, 
nous pourrons écrire : 

^+X^( F *- 7 *)+^= o (i) 

et nous aurons n équations pareilles en faisant i— 1, 2, . . . . , n. 

C'est en transformant le système (1) que Fourier est arrivé aux équations 
qui nous ont occupés dans les §§ précédents. Pour cela il passe à la limite, de fapon 
à passer des équations aux différences finies aux équations différentielles et en 
tenant compte de l'isotropie du corps. Mais si l'étude de ces équations différen- 
tielles nous conduit à une de ces difficultés qui tiennent à la considération de 
l'infini, nous ne devons pas oublier que cette difficulté est factice, puisque au 
point de vue purement physique, ces équations différentielles ne sont là que 
pour remplacer des équations aux différences finies qui en diffèrent très peu et 
pour lesquelles cette difficulté n'existe pas. Il y a donc intérêt à étudier le 
système (1) en lui-même. 

Cette étude ne présente aucune difficulté puisqu'il s'agit d'un système 
d'équations linéaires à coefficients constants. 

Posons donc : F< = U^- at , 

les Ut et a étant des constantes ; les équations (1) deviendront : 

aU^SMUt- U^+CiUi. (2) 

En éliminant les n constantes Ut entre ces n équations (2), on arrive à une équa- 
tion de degré n en a que j'écrirai : 

.F(a)=0. (3) 

Soient a lt a 2 , a n les n racines de l'équation (3); considérons une de ces 

racines que j'appelle a p . Quand dans les équations (2) on fera a = a ï >, ces 
équations deviendront compatibles et on en tirera : 
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L'intégrale générale des équations (1) sera alors : 

V 1 = A 1 e-«> t U n + Ajr«U a + • • • • + A n e~^U nl , 
V, = A 1 e-*><U n + A/r-*U n + ....+ A n e~^U nZ , 



V n =A 1 e-» t U ln + Ajr-*U»+ ....+ A n e~^U nn , 

Ai, A % , A n étant n constantes arbitraires d'intégration. 

On voit déjà par là que la véritable solution est bien de la forme à laquelle 
nous avions été conduits dans les §§2 et 3. Mais la forme symétrique des équa- 
tions (2) va nous en apprendre d'avantage. La quantité de chaleur envoyée par 
M t à M k doit être égale à celle que repoit M h de M it on a donc : 

Envisageons donc la forme quadratique positive : 

<W, u, , u n ) = xo ik (v i -v k y + xc t vî 

les équations (2) pourront s'écrire : 

Si nous regardons TJ X , JJ % , . . . . , U n comme les coordonnées d'un point dans 
l'espace à n dimensions, l'équation 

représente un ellipsoïde puisque la forme <E> est positive. 

Gomment devrait-on. procéder pour trouver les axes de cet ellipsoïde. Il 
faudrait précisément résoudre les équations (2 bis), éliminer les Ut entre ces 
équations, ce qui donnerait l'équation (3). 

Il résulte de là que l'équation (3) a toutes ses racines réelles et positives. 
Ce résultat est connue de tous sous cette forme pour les ellipsoïdes dans l'espace 
ordinaire à 3 dimensions. Il est vrai encore dans le cas qui nous occupe, comme 
nous l'apprend la théorie des formes quadratiques. 

En efifet la forme <î> peut être décomposée en une somme de n carrés et nous 
écrirons cette décomposition sous la forme suivante : 

3> = atf>\ + orfj + + <v$ 

où : %= U pl Ui + U^U t + +U pn U n , 

JJpi, Upn, etc. étant des constantes. 
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Soit maintenant : 

La théorie des formes quadratiques nous apprend que l'on peut choisir la décom- 
position de <E> de telle sorte que : 

e = tf + tf + .... + # 

ce qui entraîne les conditions : 

u P \ + 03 +.... + & P \ =1, (4) 

U p iU ql +U pg U gZ + ....+ U pn U qn = 0. (p>q) (5) 

Si donc on écrit les n équations simultanées : 

fc, = l, 4> ? =0 (p>q) (6) 

ces équations admettront pour solution : 

Vi=Upi, U,= U pi ,....U n -U pn . (7) 

Or ces équations (6) entraînent les suivantes : 

1 d<3? ___ dty p _ dty p -. J_ d® _ chpp dty p 

Y dÏÏt — a ^ p ~dU t —^dUt' Ui ~ ~2 dÏÏ i ~ $•» dût ~ dÛl 

d'où: ±_ É?t — n 

2 dU t — a p U <- 

Nous retrouvons les équations (2 bis) . Les valeurs (7) des Z7et la valeur a p de a 
nous représentent donc une solution de ces équations (2 bis). 

Quelle est maintenant la signification de ces nombres a x , a 2 , . . . . , a„. Sup- 
posons que ces nombres qui sont tous réels et positifs soient rangés par ordre de 
grandeur croissante. 

Il est clair que a x sera le minimum du rapport : 

$__ <*!<?>! + <*2<?>l + + «„«?>» 

© ~ tf + $ + .... + & ' 

Si maintenant on suppose que les n variables U ne soient plus arbitraires, 
mais soient liées par la relation 

4>i = (8) 
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on verra que a a sera le minimum du rapport : 



0- #+. 


• • • + tf 


Si les U sont liées par les deux relations 




#1 = fts = 


= 0, 


a 3 sera le minimum du rapport : 




* __ <*3<PÎ + • • 


• • + CCnfyl 



(8 bis) 



© $+.... + # 

et ainsi de suite. 

L'analogie avec l'analyse du §2 est évidente; il suffit pour retrouver cette 

analyse de passer à la limite comme l'a fait Fourier. 

Les équations (1) sont analogues aux équations de Fourier : 

Ç = «*A7, <^ + hV=0. 
dt an 

Les équations (2) et (2 bis) sont analogues aux équations qui définissent les 
fonctions U et qui s'écrivent : 

AU+hU=0, ^+hU=0. 

Les nombres a sont analogues aux nombres h. 

La forme <ï> est analogue à l'intégrale que nous avons appelée B dans le §2 
et la forme © à l'intégrale que nous avons appelée A. 

L'équation (4) est analogue à l'équation : 



fTJidv = 



et l'équation (5) (qui exprime que les axes de notre ellipsoïde sont rectangulaires) 
à l'équation: fu n D& = 0. {n> <P ) 

Il est inutile dé pousser plus loin cette comparaison, on comprend suffisamment 
la parfaite identité des raisonnements, bien que ceux-ci soient parfaitement 
rigoureux dans le cas du présent paragraphe, où l'infini n'intervient pas, et qu'ils 
soient au contraire sujets à de graves objections dans le cas du §2. 
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Ce n'est pas seulement dans l'étude du problème de Fourier qu'on est con- 
duit à ces considérations ; on obtiendrait des résultats tout à fait analogues en 
envisageant au même point de vue les autres problèmes de Physique Mathé- 
matique. 

Dans tous ces problèmes on a à intégrer des équations linéaires aux dérivées 
partielles. Ces équations ont partout la même origine. Les lois du phénomène 
véritable sont exprimées par des équations linéaires aux différentielles ordinaires, 
où la seule variable indépendante est le temps et où les inconnues sont en très 
grand-nombre; chacune de ces inconnues en effet, représente la valeur d'une certaine 
quantité relative à l'une des molécules du corps. Le nombre de ces inconnues est 
donc Jen, n étant le nombre des molécules du corps, et h le nombre des quantités 
relatives- à chaque molécule. C'est par un véritable passage de la limite qu'on 
passe ensuite de l'hypothèse moléculaire à celle de la matière continue et des 
équations différentielles ordinaires aux équations aux dérivées partielles. 

Si donc on revient momentanément à l'hypothèse moléculaire, on n'a plus 
affaire qu'à des équations linéaires ordinaires à coefficients constants, et la seule 
difficulté provient du très grand nombre de ces équations. Mais il y a plus ; ces 
équations présenteront presque toujours la symétrie que nous avons observée 
dans les équations (1) et on sera encore conduit â envisager une forme quad- 
ratique et tout sera ramené à la décomposition de cette forme en carrés. 

Je n'en donnerai qu'un exemple ; j'envisagerai les équations de l'élasticité. 
Soient x, y, z les coordonnées d'une molécule quelconque, dans l'état d'équilibre 
lorsque les forces extérieures appliquées au corps sont nulles; soient x-\-u, 
y -f v, z + w les coordonnées de cette même molécule lorsque le corps élastique 
est déformé sous l'action de forces extérieures ; soit <I> la fonction des forces rela- 
tive aux forces élastiques ; soient X, Y, Z les trois composantes de la force 
extérieure appliquée à la molécule considérée. Les équations d'équilibre s'écri- 
ront alors : $*> $*> ^j, 

~ — X, -5=- = Y, ~ = Z. (9) 

du av dw v 

Comme u, v, w sont très-petits, nous pouvons développer <ï> suivant les puis- 
sances de ces quantités et négliger les puissances d'ordre supérieure à 2. Les 
termes du 1 er degré doivent être nuls, puisque l'équilibre normal est atteint pour : 

« = v = w= 0. 

Je puis supposer que le terme tout connu est également nul ; puisque <E> n'est 
déterminé qu'à une constante près. 
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En résumé <I> sera une forme quadratique par rapport aux u, v, w, et cette 
forme sera positive parce que l'équilibre normal doit être stable. 

Les équations (9) seront donc linéaires en u, v, w. Le nombre de ces équa- 
tions est le même que celui des inconnues u, v, w; il est égal à w, si le nombre 

des molécules est — . Afin de reprendre les mêmes notations que tout à l'heure, 

o 

nous appellerons les n variables, Ci, U z , . . . . , TJ n ; alors si les trois coordonnées 
d'une molécule que nous appelions tout à l'heure w, v, w, s'appellent main- 
tenant U p , U p + i, U p+2 , nous appellerons de même X p , X p+1 et X p+2 les com- 
posantes de la force extérieure appliquée à cette molécule, composantes que nous 
appelions tout à l'heure X, Y, Z et les équations (9) deviendront : 

-^j—X p . (>=1,2 , n) (9 bis) 

Nous décomposerons la forme <E> en carrés comme nous l'avons fait tout à l'heure 
et nous retrouverons les formules : 

* = «i$ + a&l + + <VÎ>L 

<&,= u pl u x + u pi u 2 +.... + u m u n , 

© = $ +<*>! +••..+*!, 
e=u? +ui + ....+ U*.- 

Les équations (9 bis) deviennent alors : 



k — n 



J^aictyjJJhp — X P (p = 1 , 2, , n) 



(9 ter) 



fc=i 



Multiplions la première de ces équations par U a , la seconde par U i2 , la 

n 6 par U in et ajoutons. En tenant compte des équations (4) et (5) il viendra: 

«0*= TT a X l + U n X 2 + + U in X n , (•= 1, 2 , n) (10) 

Multiplions maintenant la première des équations (10) par — — , la seconde par 

U U 

— ^ . . . . , la n par — ^ ; en vertu des équations d'orthogonalité (4) et (5) ou 

plutôt des équations 

Cj p + U 2p +.... + U„p = 1 , 
UipUiq + TJ %V U H 4- .... + U np U m =0, 
38 
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qui comme on le sait leur sont équivalentes ; il viendra : 

u. = Ç±di + -^0 2 + . . . . + -^ e n 

a x a 2 2 ' a„ 

où l'on a posé pour abréger : 

Qi = ^n^i + U n X % +....+ U in X n . 

Les équations (9) sont donc ainsi résolues. 

Il est clair qu'une pareille solution ne peut être que théorique, comme 
l'était déjà la solution du problème de Fourier par l'intégration des équations (l). 
Le nombre immense des équations (1) comme celui des équations (9) s'opposerait 
absolument aux calculs. Mais cette solution purement théorique peut mettre 
sur la voie de la solution véritable. 

Passons à la limite et abandonnons l'hypothèse moléculaire pour celle de 
la matière continue. Nos équations (l) ou (9) deviendront des équations aux 
dérivées partielles ; nos formes quadratiques <î> et deviendront des intégrales 
analogues à celles que nous avons appelées A et B dans le §2. 

Notre procédé pourra s'appliquer sans autre changement ; au lieu de décom- 
poser les formes <ï> et @ en carrés, nous aurons à chercher les minima successifs 
de leur rapport ou plutôt ceux du rapport des intégrales A et B qui les rempla- 
cent. On passera ainsi d'une analyse analogue à celle de ce paragraphe à une 
analyse tout à fait semblable à celle du §2. 

§6. — Existence des Fonctions U n . 

L'existence des fonctions U„ peut être maintenant regardée comme démon- 
trée au moins au point de vue physique. La fonction Z7i sera une fonction qui 
devra prendre les valeurs 

^11 > ^12) • • • • t t/ln 

aux différents points occupés par les molécules 

M u M 2 , , M n . 

La fonction TT % devra prendre en ces mêmes points les valeurs 

Un, t/ 22 , . . . . , U 2n 
et ainsi de suite. 
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Nous savons de plus que les valeurs de la fonction TJ X par exemple devront 
satisfaire aux équations (2 bis) du § précédent que j'écrirai 

1 d ® _ tt (i\ 

-y m~ aiUi - (1) 

La fonction TJ X n'est définie ainsi, il est vrai que pour n points de l'espace, 
à savoir les n points occupés par nos n molécules ; mais comme ces molécules 
sont très nombreuses et très rapprochées les unes des autres, on pourra calculer 
par interpolation la fonction TJ X pour tous les autres points intérieurs au corps. 

On pourra à la vérité trouver ainsi deux valeurs différentes pour la fonction 
TJ X si l'on adopte deux règles d'interpolation différentes ; mais les différences 
seront du même ordre de grandeur que la distance qui sépare deux molécules et 
par conséquent négligeables au point de vue physique. 

La fonction U x ainsi définie satisfera approximativement aux équations de 

Fourier: jtt 

AÏÏ 1 + k l U l =0 1 ™l + A0i=O (2) 

que l'on obtient en partant des équations (2 bis) et en passant à la limite. 
L'erreur commise en remplaçant les équatious (1) par les équations (2) sera du 
même ordre de grandeur que la distance qui sépare deux molécules. 

Il y a à cela toutefois une condition, c'est que les dérivées de la fonction TJ X 
soient finies ; on n'aurait plus le droit de passer des équations (l) aux équations 
(2) étaient du même ordre de grandeur que l'inverse de la distance qui sépare 
deux molécules. 

Il nous resterait donc à établir que ces dérivées sont bien finies; c'est-à-dire: 

1°. Que la différence U u — U lk est du même ordre de grandeur que la 
distance des molécules M t et M K . 

2°. Plus généralement, soit M a une molécule quelconque de coordonnées 
x, y, 2, soient M p , M y , . . . . , M K un certain nombre de molécules très voisines 
de M a et dont les coordonnées soient respectivement, 

* + &» y + vp, * + &; » + £ y , y + ri y , z + £ y ; — /» + &, y + nx, s + &. 

Soit P(Jz, y, Ç) un polynôme quelconque de degré inférieur àmenf, y, Ç. 
Soient enfin a, (3, y, . - . . , /l un certain nombre de coefficients relatifs aux 
diverses molécules M a , M ? , M K ; et supposons que ces coefficients satisfas- 
sent à la condition suivante : 

<*P(£«, >?., Q + (3P&, v„ Q + yP&, n y , Q + • . • . + *P&, n„ Q = o 
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et cela quel que soit le polynôme P pourvu que son degré soit inférieur à m 
(j'écris par symétrie P(£ a , y; a , Ç a ) au lieu de P (0,0,0)). 

Si ces conditions sont remplies, nous aurions à établir que : 

est du même ordre de grandeur que les puissances m es des quantités £, n, £. 

Cela ne serait sans doute pas impossible ; je n'aurais en effet pour établir ces 
divers points qu'à traduire dans le langage de l'hypothèse moléculaire l'analyse 
du §4. 

C'est ce que je me réserve de faire dans un mémoire ultérieur qui pourra 
être regardé comme la suite de celui-ci. 

Je pourrai dire alors que les conclusions des §§2, 3 et 4 sont démontrées d'une 
fapon rigoureuse au point de vue physique. Peut être même est-il permis 
d'espérer que, par une sorte de passage à la limite, on pourra fonder sur ces 
principes une démonstration rigoureuse même au point de vue analytique. 
Paris, le 19 Mars, 1889. 



